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Abstract 

Fermilab Experiment 866 took single muon data from 800 GeV/c {y/s = 38.8 
GeV) p-Cu and p-Be interactions in an attempt to extract the inclusive nuclear open 
charm/anti-charm {D/D) differential cross sections as a function of pr- The muons 
were decay products from semi-leptonic decays of open charm mesons as well as decays 
from lighter non-charmed mesons (tt's and K's). Data were taken simultaneously from 
two interaction regions; one of two thin nuclear targets and a copper beam dump 92 
inches downstream. The open decay length for hadrons produced in the targets 
increased the contribution to the muon spectrum from light hadron decays, relative 
to those from the dump. Production cross sections for light hadrons from previous 
experiments were used in conjunction with parameterized open charm cross sections 
to produce total Monte Carlo single muon spectra that were subsequently fit to the 
data. 

The sensitivity of this measurement covered an open charm hadron pt range 
of approximately 2 to 7 GeV/c, center-of-mass rapidity, Ucm, between and 2, and 
xf between 0.2 and 0.8. Previous experimental results for p-p or p-A open charm 
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production at comparable energy was limited to -\/5 GeV/c. Three functions describ- 
ing the shape of the open charm/anti-charm cross sections were fit to the data; an 
exponential, Aiexp{^BpT), and two polynomials, ^ ^ and A2 ^^J^^,[J^^'"'°^n . 

The first polynomial was fit with the parameter n as a free parameter, and constant 
with three integer values, 4, 5 and 6. The second was fit with n held fixed at the 
constant integer values only. The best results were with the first polynomial with 
n around 6. All three parameterizations resulted in good fits. Extrapolation of the 
cross sections to small Pt shows good agreement with previous experiments. The 
power a of the nuclear dependency A°'^'^'> was calculated as a function of pr- The 
result indicates that a is transverse-momentum dependent, albeit within large errors. 
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Chapter 1 
Introduction 



Quantum Chromo-Dynamics (QCD), and the Standard Model (SM) as a whole, have 
been remarkably successful in describing the nature of particles and their interactions. 
They are widely accepted by the community. The Standard model is composed of six 
quarks (referred to as flavors), six leptons and four force carriers, neglecting gravity. 
Quarks and leptons come in three generations. A quark generation consists of two 
quarks, one having charge + 2/3 and one having charge —1/3. The quarks also 
come in three colors, red, green and blue. Each generation of leptons consists of a 
pair with one neutral and the other having unit charge. Quarks and leptons have half- 
integer spin and couple (depending on the interaction) to the spin 1 force carriers, 
gluons, the charged and neutral weak bosons {W~^ , W~ and Z^) and photons. Table 
11.11 shows the quarks, leptons and gauge bosons used in the Standard Model, taken 
from the Particle Data Group jjlj. 

Hadrons composed of three quarks, such as the proton {u, u, d) and neutron 
(m, d, d) are referred to as baryons, and hadrons made of one quark and one anti-quark 
are called mesons. The (u, u, d) quarks in the proton are referred to as valence quarks. 
Hadrons have virtual qq pairs as well, which are referred to as sea quarks. Flavor 
is conserved in strong and electromagnetic interactions, so net 'new' flavor in strong 
or electromagnetic interactions must be 0. Charm, bottom and top quarks are often 
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referred to as heavy quarks. Mesons containing charm come in two varieties, those 
having one charm or anti-charm quark, referred to as open charm, such as the 
which has one charm and one anti-down quark [c dj , and those containing one charm 
and one anti-charm quark, referred to as hidden charm, such as the J / \1/ (cc). 

Study of the production of hadrons containing charm or heavier quarks con- 
tributes to the understanding of the theory of Quantum Chromo-Dynamics (QCD). 
QCD uses factorization theory to describe the production and hadronization of heavy 
quarks into hadrons seen in the lab. Figure 11.11 shows a parton-parton (partons are 
quarks or gluons) interaction in a colhsion of two baryons, A and B. The interac- 
tion produces one charm and one anti-charm quark, c and c. The charm quark then 
hadronizes into a charm meson, shown as in the Figure. 




Figure 1.1: Representation of qq cc from the interaction of a quark q in hadron A 
with momentum fraction Xq and a sea anti-quark q from hadron B having momentum 
fraction Xg. The charm quark subsequently hadronizes to a meson at the process 
labelled F^. 
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A shorthand method used to describe what is shown in Figure 11.11 is: 

A + B^D+ + D-+X 

where the charm quark hadronized into a D~^, the anti-charm quark hadronized into 
a D~ and X is used to denote 'anything else', which must include two baryons. 

All mesons are unstable and decay. The short hand notation for one of the 
ways a meson may decay, called a hadronic mode, is: 

D+ K- 7r+ 7r+ 

An example of a decay mode referred to as a semi-leptonic mode is 

Experiments studying the production of open charm usually 'tag' events by looking 
for the decay products from various decay modes. For the hadronic mode above, the 
event would be tagged if a K'^ and 2 vr"'" tracked back to a common vertex, suggesting 
they were the decay products from the . The invariant mass is then calculated, 
and compared to the mass of the . If the event passes all cuts used to eliminate 
bad events, it is used in the analysis of the production characteristics of the for 
their experiment. 

This analysis used data from 800 GeV/c p-Cu and p-Be interactions to de- 
termine the production of hadrons containing open charm as a function of the hadron 
Pt- The data consisted of single muon events. The E866 spectrometer was designed 
to study di-muon events, where the two muons tracked back to a decay vertex. Single 
muon events are events where a hadron decays semi-leptonically to one muon, plus 
anything else, as shown above. The E866 spectrometer could not track any remain- 
ing secondary particles from the decay, so no decay vertex was available. The single 
muon data was taken from two production regions simultaneously, one a thin target 
of either copper or beryllium, the other a solid copper beam dump. The targets 
had an open decay length of 92 inches before hadrons (and remaining proton beam) 
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interacted in the dump. Hadrons containing open charm decay roughly IC^ times 
sooner than hght hadrons, so the ratio of muons from open charm to hght hadrons 
is significantly enhanced in the data taken from the dump, relative to the data taken 
from the targets. Use of the difference between muon spectra from the two production 
regions allowed extraction of the inclusive open charm/anti-charm differential cross 
sections as a function of the hadron transverse momentum between 2.25 and 10.0 
GeV/c. Mesons containing open charm are shown in Table IT!^ 

The term inclusive (versus exclusive) is used to classify the production of open 
charm by the final state of the reaction 

p + N ^D + X 

where p is a proton interacting with a nucleon in the material, D is a hadron 
containing open charm or open anti-charm and X is anything else.^ This analysis 
could not identify the parent hadron, so any hadron produced that decayed to an 
open charm or anti-charm hadron is included in the measured cross sections as if it 
were produced as an open charm or anti-charm hadron, such as 

p + N ^ B + X where B D + n 

Exclusive cross sections measure the strength of a reaction to a specific set of particles 
such as: 

p + N^TT^ + TT~+p + N 

Additional information can be derived as well, such as the production depen- 
dency on A."^ Prior experimental results from meson interactions have shown an 

"'^To conserve baryon number, X must be composed of at least two baryons and the other charm 
quark must be included in either the baryons or as another meson. 

is the atomic weight of a material used as a target to produce hadrons in interactions. The 
scaling of the production to the atomic weight is commonly referred to as the nuclear dependency 
which is discussed in Chapter^for the production of light hadrons, as well as in Chapter regarding 
the values for open charm production as determined from this analysis. 
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enhancement of the production of open charm hadrons containing one of the valence 
quarks of the incident meson, called the leading particle effect. This enhancement has 
not been observed in proton interactions. If this effect were present in proton-nucleon 
interactions, it would be seen as an enhacement of hadrons relative to anti-hadrons. 
The hadron/anti-hadron ratio determined by this analysis is presented as well. 



1.1 Heavy Hadron Production 

The theoretical description of heavy hadron production in pQCD (pcrturbative Quan- 
tum Chromodynamics) is done in two parts, referred to as factorization; production 
of the heavy quarks using partonic cross sections, and the process of hadronization, 
where the bare quarks are transformed into hadrons seen in the lab. 

A cross section, cr, is used to measure the effectiveness of an interaction such 

as 

a + 6 — > c + d 

where a and b are the interacting (incident) particles and c and d are particles pro- 
duced in the interaction which may be different from cither a or b. The cross section 
has units of area, typically given as barns (1 barn = 10~^^ cm^). 

Partonic cross sections in QCD are modeled from Deep Inelastic Scattering 
cross sections. The cross sections for deep inelastic scattering on unpolarized nucleons 
can be written, generically, in terms of structure functions: 



+ (1.1) 

where i = NC,CC is for neutral or charged current scattering and Fl, and Fl 
are structure functions. Q is the four-momentum transferred in the interaction, and 
a (called the fine structure constant) is defined as a — e/tic. (The units for are 
(GeV^ cm^).) In the quark-parton model, contributions to the structure functions 



(Pa' _ ina' 
dx dy xyQ"^ 
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Table 1.1: Quarks, leptons and vector bosons used in the Standard Model. Masses 
(MeV c~^) are given in parentheses pQ. 



Quarks and Leptons 


Type 


charge 


1'^* Generation 2^'^ Generation 3'"'^ Generation 


up- type 
quarks 


2 
+ 3 


u c t 
up charm top 
(1.5 - 4.0) (1150 - 1350) (174000) 


down-type 
quarks 


1 

~3 


d s b 
down strange bottom 
(4 - 8) (80 - 130) (4100 - 4400) 


neutral 
leptons 





I/e Vr 

e neutrino /i neutrino r neutrino 
(< 3x10-6) (<0.19) (< 18.2 ) 


charged 
leptons 


+1 


e /i r 
electron muon tau 
(0.5) (105.7) (1777) 


Vector Bosons (Force Carriers) 


Type 


charge 


Boson 


mass 




neutral 
weak 





Z 


(9120) 


charged 
weak 


+ 1 
- 1 


w- 


(8040) 


electro- 
magnetic 





7 

photon 


(< 6 X 10-23) 


strong 





9 

gluon 


(0) 
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Tabic 1.2: Mesons containing open charm. Mass (MeV/c^) is rounded off to tlie 
nearest MeV/c^. 



Pscudoscalar Mesons (Spin 0) 




D+ 




Dt 


D- 


If 


D- 


Quarks 


c, d 


c, u 


c, s 


c, d 


c, u 


c, s 


Mass 


1869 


1865 


1969 


1869 


1865 


1969 


Vector Mesons (Spin 1) 






D* + 


D*^ 


D*- 


D*' 


Quarks 


c, d 


c, u 


c, d 


c, u 


Mass 


2010 


2007 


2010 


2007 



7 



are expressed in terms of quark distribution functions of the proton. The quark distri- 
bution function, q{x,Q'^), is the number of quarks, or anti-quarks, of the designated 
flavor {q = u,u,d,d, s, ■ ■ ■) that carry a momentum fraction between x and x + dx 
of the protons momentum where the protons momentum is large. For the charged 
current reaction e~p —>■ vX 

=2x{u + d + s + c+---) 

One prediction of the quark-parton model is that the structure functions scale 
in the Bjorken limit that and z/ — oo with x fixed. Scaling implies 

F\x,Q^)^F\x) 

QCD uses scale dependent parton distribution functions (PDFs), /(x, /x^), to describe 
the process above. Here, f = goiq and fi is typically the scale of the probe, Q. At a 
given X, these correspond to the density of the partons in the proton integrated over 
transverse momentum kt up to fi/c. The fine structure constant a used in Equation 
II. H is redefined to be scale dependent Q;s(/i^) = ^ where gs is the SUc(3) coupling 
constant. Parton distributions have been measured by many experiments. Figure IT!^ 
shows the distributions x times the unpolarized parton distribution functions f{x) at 
a scale /i^ = 10 GeV^/c^^. 

Production of heavy quarks in hadron-hadron interactions at leading order 
(LO) in perturbation theory is the result of two parton-parton interactions 

qq~^QQ and gg^QQ 

where qq (often referred to as quark-quark annihilation) represents a quark from one 
hadron interacting with an anti-quark from the other hadron, g g is used to represent 
the interaction between a gluon in one hadron interacting with a gluon in the other. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



X 

Figure 1.2: Distributions of x times the unpolarized parton distributions f{x) (where 
/ = Uy, d^,u, d, s, c, g) using the MRST2001 parameterizations |21IS1 (with uncertain- 
ties in Uy, dv and g) at a scale /x^ = 10 (GeV^/c^). Figure is taken from p. 

and Q Q is the produced heavy quark and anti-quark (they must be produced as a 
pair). In p-p or p-A colhsions, qq interactions are between a valence quark in one 
hadron, and a sea anti-quark in the other. 

Feynman diagrams are used to describe the process as well as to calculate 
the amplitudes of the interaction. The Feynman diagrams of the leading order (LO) 
processes for the production of charm are shown in Figure 11.31 The total partonic 
cross section from hadro-production is proportional to the sum of all the combinations 
of the quarks in one hadron with anti-quarks in the other (and conversely, the sea 
anti-quarks in the first with quarks in the other) plus the contribution from gluons 
in both. 

Hadronization is usually described in terms of fragmentation functions, -F(x, s). 
The Particle Data Group PP define the fragmentation function for a hadron of type 
h at cm. energy a/s represented as a sum of contributions from the different parton 
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Figure 1.3: The Leading Order Feynman diagrams qq ^ Q Q (1) and g g ^ QQ (2). 
types i — u,d,s,u,d,s, - ■ ■ , g' as 



where D!^ are the parton fragmentation functions analogous to the parton distribution 
functions above, x — 2 E^/ s < 1, z — x^/xi and Cj are coefficient functions of the 
partons i. 

At lowest order in a^, the coefficient function Cg for gluons is 0, but for quarks, 
Ci = gi{s) 6 {1 — z). At higher orders of as the coefficient and parton fragmentation 
functions are factorization scheme dependent. Measured fragmentation functions need 
to be parameterized at some initial scale, to, typically 2 GeV^ for light quarks and 
gluons. A general parameterization is: 



where the normalization N ^ a, (3 and 7 are usually dependent on the scale to as well 
as the type of parton, p and hadron h. 






1 
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There have been numerous studies of the production of open charm using 
nuclear targets. The majority of fixed target studies have been done using charged 
meson beams on nuclear targets. A reasonably complete compilation of experimental 
results is found in Here, and in most literature, the term open charm production 
is used as a 'catch all' phrase for all hadrons containing one charm or anti-charm 
quark. Results are usually presented for the total open charm production, or one or 
more sub-categories of these hadrons, such as hadrons containing one charm quark 
{D) or one anti-charm quark (D) or those with charge ±1 or neutral. 

1.2 Results From Prior Experiments 

Discussion concerning results from fixed target p-p and p-A open charm production 
often cite publications from four experiments; Fermilab Experiment 743 (the LEBC- 
MPS Collaboration) Fermilab Experiment 653 |2|, Fermilab Experiment 789 fH] 
and Fermilab Experiment 769 ■ Figure 11.41 shows the total open charm differ- 
ential cross section from 800 GeV/c p-p interactions as measured by the LEBC-MPS 
Collaboration jB]. 

Fermilab Experiment 789 measured neutral open charm production for 800 
GeV/c protons incident on beryllium and gold targets ^8J. Their results are shown in 
Figure IT31 The Fermilab E769 Collaboration studied open charm production using 
secondary 250 GeV/c vr^, and p beams on a multifoil target of beryllium, copper, 
aluminum and tungsten j^lE]. The results from are presented in Figure EH 

The production of open charm using charged meson beams shows significant 
differences from the production of open charm using proton beams on nuclear targets. 
The differences are seen in the ratio of charged to neutral charm production and the 
ratio of hadron to anti-hadron production. The second is referred to as the leading 
particle effect. Results to date from either meson or proton induced production reveal 
that open charm meson production is shaped more or less the same as the next-to- 
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leading order (NLO) predictions for the production of the quarks themselves (see 
Figure II .bp . 

Theory suggests that the ratio of charged to neutral open charm production 
should be approximately 0.32 [SJ and results from open charm production us- 
ing charged meson beams shows rough agreement. The results from fixed target 
p-p and p-A experiments, until recently, have shown this ratio to be more or less 
unity. A recent result published by the HERA-B Collaboration reports the ratio 
a{D+) / (t{D^) = 0.54±0.11±0.13 from 920 GeV/cp-^ induced production y^. The 
STAR Collaboration reported the same ratio as 0.40 ± 0.09 ± 0.13 using data from 
^/sNN = 200 GeV p-p and d-Au collisions [T^. The leading particle effect reported 
by meson-nucleon experiments is contrary to predictions as well. Proton-nucleon ex- 
periments have not seen this effect. It is generally thought that some momentum is 
lost during the hadronization process. Results from a variety of experiments, such as 
E769 above, do not show this effect. 
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Figure 1.4: Total inclusive open charm differential cross section da{D + D) / dp^ 
(yub c^ GeV~^) measured by the LEBMC-MPS Collaboration (Fermilab Exper- 
iment 743) [6] • Solid curve shows the results of a fit to the empirical form 
(1 — exp {—ap^). 
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(GeV/c) 



Figure 1.5: Differential cross section per nucleon, da{D^) / dpT + da{D ) / dpT versus 
Pt- Uncertainties shown are statistical only and do not indue an additional systematic 
uncertainty of 12.8 percent. Figure is from Fermilab Experiment 789[H]. 
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Figure 1.6: Measured D meson (£)+, D-, D , D+andD;) da / dpi 
(yubc^ nucleon"^ GeV~^) {xp > 0) for production induced by tt, K and p beams and 
NLO QCD predictions for charm quarks [TT] (tt and p beams). In addition to the 
statistical errors shown, there are overall normalization errors of about 6%, 6% and 
9% for TT, K and p results respectively. Figure is taken from jH]. 
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Chapter 2 



Apparatus 

2.1 General Description 

Fermilab Experiment 866 (E866/NuSeA) was designed for collecting and analyzing 
di-muon events from 800 GeV/c protons incident on various nuclear targets. E866 
was a continuation of several Fermilab Experiments including 772 and 789, where 
the detector had several major improvements to increase the accuracy of the mea- 
sured trajectory of the muons as well as increased data taking capabilities through 
improvements in the data aquisition system and trigger configurations. The exper- 
iment has resulted in five previous Doctoral theses JH] - 120] and several published 
results Ell - Ell- 
Figure shows the FNAL E866/NuSeA spectrometer for the original config- 
uration of the experiment. There were modifications to the original configuration to 
take data for this analysis that will be pointed out and explained as necessary. 

The E866 spectrometer was located in the East hall of the Meson Experimental 
Area, and used the 800 GeV/c proton beam extracted from the accelerator ring for 
a period of approximately twenty seconds out of each minute. Each spill contained 
"buckets" with a 19 ns time structure based on the radio accelerator RF of 53 MHz. 
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Figure 2.1: FNAL E866/NuSeA Spectrometer. For this analysis the spectrometer had 
no cryogenic target system, but a target wheel instead, missing is the dump shown in 
figure 1221 The ring-imaging cherenkov counter was inactive, and the muon detectors 
are referred to as station 4 in the text. 



Beam monitoring was accomplished with various detectors for size, position and in- 
tensity. The beam size and position was measured using RF cavities and segmented 
wire ion chambers (SWIG). The beam was last monitored approximately 70 inches 
upstream of the targets. This was done using a movable SWIG having a wire spacing 
of 2 mm horizontally and 0.5 mm vertically. There were several intensity monitors, 
but these were unreliable at the low intensity (approximately 10^° protons per spill) 
used for the single muon data. 



The data was taken with a rotating target wheel instead of the cryogenic target 
system shown in Figure 12.11 The experiment designated the center of the opening of 
the upstream face of spectrometer magnet SM12 as the origin. The target wheel was 



2.2 Beam Monitors 



2.3 Targets 
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located aX z = —24 inches. The wheel had three used positions; one empty target 
frame, one beryllium target of thickness 2.036 inches (referred to as the beryllium 
target), and another copper target of thickness 1.004 inches (referred to as the copper 
target). The proton interaction lengths of these targets, calculated with pcu = 8.96 
gm cm~^ and Xi^cu = 134.9 gm cm~^ for copper and pse = 1-848 gm cm~^ and Xi^Be = 
75.2 gm cm~^ for beryllium, were 0.127 and 0.169 interaction lengths respectively. 
The wheel was rotated during the 40 second period between spills according to two 
rotation schedules. 

2.4 Spectrometer 

The spectrometer had three dipole magnets with fields in the x direction as defined in 
Figure 123 For this analysis data was taken with currents in magnets SM12 and SM3 
only. The currents were configured both parallel and anti-parallel during data taking, 
although results are only shown for the parallel configuration. The magnet currents 
were set to 1420 amps in SM12 and 4200 amps in SM3. Only the results from the 
parallel magnetic field configuration are given in this analysis, because the analysis 
lost the computer during the analysis of the opposite polarity data. The magnetic 
configuration would only allow for a separate set of data to fit, since combining the 
two results was not possible. The acceptances for both magnetic field configurations 
was very similar, so no new information was lost. 

Particles created in the targets, and any remaining proton beam, entered mag- 
net SM12. The configuration of this magnet when the data for this analysis was 
taken is shown in detail in Figure 17!^ The magnet was 567 inches long and provided 
a momentum deflection of 7 GeV/c when operated at its maximum current of 4000 
amps and provided a momentum deflection of approximately 2.4 GeV/c for this data. 
The copper beam dump was 168 inches long beginning 68 inches inside the volume 
of the magnet. The dump spanned the entire volume in x, and was approximately 8 
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inches tall for the first 80 inches, and 10 inches tall thereafter. The face of the dump 
had a rectangular hole 12 inches deep by 2 inches square to reduce back-scattered 
particles. The dump provided a second target having 26 proton interaction lengths, 
or approximately 220 radiation lengths. 

Behind the dump was the hadron absorber that filled the entire volume in 
X and y consisting of 24 inches of copper, 4 layers of carbon 27 inches thick each 
and two 36 inch layers of borated polyethylene. The absorber wall, having over 
thirteen interaction or sixty radiation lengths, effectively allowed only muons to pass 
downstream. The remainder of the inside of the magnet was filled with a helium bag. 

Between magnets SM12 and SMS is the first of three similar tracking stations. 
Stations 1, 2 and 3 had multiple pairs of drift chambers and one or more layers of 
hodoscopes. Station 1 had three pairs of drift chambers, Ul, Ul', Yl, Yl' and VI 
and VI'. Planes U and U' had sense wires oriented +14° [tan{6) = 0.25] in the x — y 
plane, Y and Y' were horizontal and V and V were oriented -14°. The primed planes 
were offset in the direction perpendicular to the sense wires by half of one cell to help 
remove ambiguities in the drift direction. Table lA.ll (page I142p gives a summary of 
the physical characteristics of the drift chambers for each station. 

Hodoscopes at each station allowed fast track evaluations or 'roads' used to 
trigger valid events to be taken to tape. Horizontally aligned planes of scintillators 
determined rough y positions, while vertically aligned planes gave rough x positions. 
Each plane was optically split to provide quadrants named up left (UL), up right (UR), 
down left (DL) and down right (DR). Each plane was designated by orientation and 
station number. Gaps created in the splitting resulted in small dead spots for muon 
tracks at or near a; = for y hodoscope planes, or y = for x hodoscope planes. 

Specific to this analysis, data were taken with the middle half of each x mea- 
suring hodoscope plane turned off to reduce trigger rates that would have been un- 
acceptabes had they been left on. This configuration was necessary to reduce the 
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event rate, even at the low intensity requested. This configuration produced a loss in 
acceptance for muons having px (/u) < 1.75GeV/c, and created two inner acceptance 
edges which further reduced the acceptable minimum muon transverse momentum as 
described in I2.5l and l3.3l Figure ESI shows an x hodoscope plane, and Table gives 
the specifications of all hodoscope planes as used for taking data for this analysis. 

Particles then entered the second magnet, SM3, which was used to determine 
the momentum of the particle in the y — z plane. At its normal current setting of 
4230 Amps, the magnet provided a transverse momentum defiection of 0.91 GeV/c. 
The direction of the defiection relative to the magnetic field direction determined the 
charge of the particle. Muons then passed Station 2, a tracking station similar to 
Station 1 except only one plane of y measuring hodoscopes, Y2, is present. 

Muons then traversed a non-operational Ring Imaging Cherenkov counter 
(RICH). This detector had been used in previous experiments for particle identifi- 
cation, but was not used for this analysis. The interior was helium filled to reduce 
multiple scattering. Downstream of the RICH detector was Station 3, a larger version 
of the prior stations. Station 3 had both x and y measuring hodoscope planes. 

Both calorimeters following Station 3 shown in Figure 12.11 were inactive, but 
had been left in place to provide additional hadron absorbing material. Shown in 
Figure 12.11 are large amounts of zinc and concrete placed between three layers of 
proportional tubes, PTYl, PTY2 and PTX in Station 4, labeled as 'Muon Detectors' 
in figure ITTl Station 4 had two hodoscope layers for measuring both x and y, called 
Y4 and X4. 

All wire planes and proportional tubes used the same gas mixture, 50% argon, 
50% ethane and a small amount of ethanol alcohol added by bubbling the gas through 
ethanol which was kept at a constant 25° F. 
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Figure 2.2: Spectrometer magnet SM12 as configured for tliis analysis sliowing tfie 
copper beam dump and fiadron absorbers in plan view (top) and elevated view (bot- 
tom). Cross indicates tfie location of the targets relative to the spectrometer magnet 
and dump. 
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Figure 2.3: View of an x plane of hodoscopes. Gaps between quadrants are exagger- 
ated. Only one phototube is shown. Shaded region represents hodoscopes having high 
voltage supplied to the attached phototubes, y hodoscopes are similar except being 
rotated 90 degrees and all scintillators have high voltage supplied to the phototubes. 
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2.5 Special X Hodoscope Setting 



The single muon analysis used a special configuration of the x measuring hodoscopes. 
The high voltage supply for the middle half of all x hodoscope layers was turned off 
(referred to as being pulled), shown in figure 17!^ This reduced the event rate for single 
muon events as well as providing a way to distinguish between events from the target 
and dump. This was a significant change from the normal operation of the apparatus. 
Single muon events required hits on the same side in all three x hodoscope layers. 
This X hodoscope configuration plus the required trigger, described in 12. 61 limited the 
transverse momentum of accepted events, because the muon must have a minimum 
absolute value for the slope in the x — z plane, |TAN^^| = [p^xl^zl- This minimum 
slope is shown in figure 123 The figure also shows how the loss in acceptance allowed 
separation of target and dump events for single muons. The minimum and maximum 
|TAN6'a;| for muons to be accepted, based on the spectrometer survey are listed in 
Table 



2.6 Event Triggers and Readout 

The FNAL E866 /NuSeA trigger has been described in detail in several references j2H] 
PU] . The data taken for the single muon analysis contained both single and dimuon 
events, requiring dimuon as well as single muon triggers. Dimuon events were used 
to calibrate the analysis routine as described in Chapter El The calibrations found 
from the di-muon events were then used to analyze the single muon events. 

Signals from the hodoscopes were first sent to LeCroy 4416 discriminators 
whose signals were reshaped and synchronized to the accelerator RF clock. All triggers 
are hit correlations from some or all hodoscope layers and are referred to as Physics 
A, Physics B and Diagnostics triggers. The single muon data was taken using 5 single 
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Figure 2.4: Drawing of change in acceptance of the spectrometer in the x — z plane, 
as well as the minimum and maximum opening angles caused by disconnecting the 
high voltage supplies to the middle half of the hodoscopes in all x hodoscope planes. 
Note the extreme shortening of the z axis, though relative distance from targets 
to hodoscope planes is preserved. Maximum angle refers to the maximum physical 
opening, while maximum angle allowed is the angle in the final cuts to insure all 
events passed through a hodoscope in all three planes. 
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and 5 dimuon physics triggers in two separate but similar configurations. The trigger 
system for the left half of the spectrometer is shown in figure lA.ll [TB] . 

The trigger system had a total propagation delay from inputs from hit ho- 
doscopes to output of the event to memory storage of approximately 20 nsec, pro- 
viding single bucket resolution. This was accomplished by comparing hodoscope hits 
in the seven layers to hardwired and software defined matrices in modules called 
Track Correlators (TCs) or Matrix Modules (MMs). The four Trigger Matrix mod- 
ules (called MUL, MUR, MDL and MDR) were lookup tables loaded onto a set of 
ECL SRAM chips. These modules used correlations in the Yl, Y2 and Y4 hodoscopes 
only. A track of interest would define a 'road' in the Y view under the deflection of 
the magnetic fields in SM12 and SMS. Hits on Y2 and Y4 were correlated against 
predicted hits on Yl and if the desired correlation existed, output for the event was 
generated. These were used only for dimuon events. Four S4XY Track Correlators, 
also used for di-muon events only, used signals from the X4 and Y4 hodoscopes only. 
The correlations required for a di-muon event to be accepted could be changed by 
software. 

The three x hodoscope layers were hardwired into an and for each side, called 
X134L and X134R. Any muon which traversed an energized hodoscope paddle in all 
three x layers on one side satisfied the XI 34 trigger. 

These were passed to the two main physics Track Correlators (along with 
the other information from the Matrix Modules and S4XY Track Correlators) called 
Physics Triggers TC and Diagnostic Triggers TC. These two main correlators deter- 
mined the validity of an event and if it was to be taken to tape, then informed the 
Master Trigger to set the busy signal and stream the event to a large memory buffer, 
and subsequently written to 8mm tape. 

The di-muon triggers were labeled Physics TcAl through 4, Diagnostic 3 and 
all four S4XY triggers. All dimuon triggers were not prescaled due to the low beam 
intensity. 
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One example of an opposite sign dimuon pair physics trigger was a coincidence 
of a track through the upper and lower left (or alternatively right) quadrants (one 
being the /i"*" and the other the and an example for a like sign dimuon event 
where one muon goes through the upper right quadrant in coincidence with a muon 
passing through the upper left quadrant, all using the four Matrix Modules MUL, 
MUR, MDL and MDR. The Diagnostics trigger required a left-right coincidence in 
five of the seven layers of hodoscopes on both sides of the spectrometer. 

The X134L/R three-fold coincidence trigger, called PhysBl, was the exclusive 
single muon trigger used for analysis. For diagnostics purposes there were 4 other 
single muon triggers, PhysB2, Diagl, Diag2 and Diag4. The diagnostic triggers were 
useful for studying edge effects of the pulled hodoscopes as well as hodoscope efficien- 
cies for two of the three X layers. 

The first trigger, PhysBl, required a hit in all three X hodoscope banks on 
either side. PhysB2 was a prescaled trigger to take 1 of 2000 events where there was 
a track through a quadrant (any one quadrant of the four) using the MUR/L MDR/L 
matrix modules, thus avoiding the large opening angle required for trigger PhysBl. 

The main single muon event trigger was PhysBl, called the X134L/R trigger 
referring to the three fold coincidence of the three X hodoscope planes on either side, 
left or right. This trigger is often referred to by the single muon analysis as the 
'trigger bit 5' trigger since that bit was set if the event had that coincidence. All 
single muon events were required to have a valid PhysBl trigger set. 

All events satisfying one or more triggers were written to tape. This was accom- 
plished using Nevis Transport electronics [21], and VME as well as CAMAC modules. 
All detector subsystems fed data onto the Transport bus, which was then transferred 
to the memory buffer. The memory buffer, which was VME based, formatted the 
event data and transferred it to 8mm tape. Since the typical spill lasted 20 seconds 
out of every minute, the memory buffering allowed higher event rates with less dead- 
time. Each event written to tape contained the values from the coincidence registers 
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from the hodoscopes and proportional tubes and readouts from the time to digital 
converters (TDC's) from signals in the drift chambers, as well as beam position, size 
and intensity, target position, voltages and information from various monitors. 

2.7 The Single Muon Data 

Table 2.1: Data taken for this analysis. All data were taken with a common trigger 
configuration. Currents in SM12 and SMS set to -1400 and -4200 amps, respectively. 
All data was taken using the same trigger configuration file, labeled jpsism2. The 
triggers for this configuration are given in the text. 



Run 


Date 


Time 


Events 
X 10^ 


IVigger 


2748 


07/10/97 


10:14 


5.778 


jpsism2 


2749 


07/10/97 


14:06 


6.047 




2750 


07/10/97 


15:41 


6.067 




2751 


07/10/97 


17:16 


5.918 




2752 


07/10/97 


19:40 


5.942 




2753 


07/10/97 


21:42 


5.750 




2754 


07/10/97 


23:30 


7.995 




2755 


07/11/97 


01:38 


6.172 




2756 


07/11/97 


03:41 


5.944 




2757 


07/11/97 


06:06 


3.150 
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Chapter 3 
Analysis 



The 10 runs on 10 tapes listed in Table ITT] fpage I^B j) contained the single muon data 
used in this analysis. The E866 collaboration analyzed data for dimuon events with 
an analysis routine first developed for the E605 experiment. The dimuon code was 
modified to analyze both single and dimuon events for this analysis. All modifications 
to the original routines are presented in detail. Analysis of the single muon data was 
performed in four stages; 

1. Separate single from dimuon events, separate target events from the dump 
events, then place the separated raw data into Data Summary Tapes, or DSTs, 
for further analysis. 

2. Iterative analyses of the dimuon DSTs were then performed to determine the 
correct values, or calibrations, needed to initialize the code. No further analysis 
of dimuon data was performed. 

3. Multiple analyses of the single muon data was performed using the single muon 

DSTs, placing the fully analyzed events into large arrays called n-tuples.^ 

^The n-tuple as well as PAW which was used to do the final analysis of the data were produced 
by CERN. 
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4. Output final data (the single muon spectra) to arrays for use in the fitting 
routines to determine the open charm cross sections. 

Figure 13.11 shows a fiow chart of the method used to analyze the single muon 
data, and may be useful for the following discussions. All analysis of the raw data 
as well as generation and analysis of Monte Carlo events was performed on a DEC 
Alpha 500 workstation. 



3.1 Initialization and Unpacking 

The first step in the analysis code was to initialize for the specific run. This initial- 
ization included setting the magnetic fields in SM12 and SMS to the current direction 
and value at the start of the run, trigger information and input spectrometer cali- 
brations. The code identified the number of the run from the data tape and used 
an internal look up table to set the magnetic field configurations whereas the trigger 
information and spectrometer calibrations needed were read from specified trigger or 
data tables. 

The code then unpacked and processed the events. Unpacking an event was a 
process that read the output DAQ information contained on the tape for that event 
and allocated the information to several large arrays which were designed to speed 
up the analysis of each track. 

3.2 Tracking 

Tracking began with searching for hit clusters in all the wire planes in Stations 2 and 
3 (See Figure Em page IT^ . This required hits in 4 of the 6 planes in each station to 
qualify as a potential hit cluster. Those clusters having hits in 2 of the 3 planes were 
classified as doublets, those having all three planes registering a hit were classified 
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Figure 3.1: Flow chart of the process used to analyze the single muon data. The data 
was histogrammed as a function of the muon pt and placed into arrays for use in the 
least-squares minimization routines described in Chapter |S] Each target resulted in 
four spectra, one each for /i^ and fi~ from the target and one each from the dump. 
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as triplets. All clusters required at least one hit in the Y view plane. Once a hst of 
doublets and triplets had been made, combinations of hits in the two stations were 
linked to produce track candidates. The linkage was a combination of hit clusters 
that loosely pointed back to the target region. 

The next tracking process was to extend the candidate tracks to hit clusters 
in the Station 1 chambers. Each candidate track was projected upstream through 
SM3 in the x-z plane to a vertical band in Station 1. Only clusters found within 
this band were considered for that track. If a cluster was found in the band, the 
entire track from Station 1 to Station 3 was refit using all 18 wire planes in the three 
stations. The z coordinate in SM3, Zsms was allowed to vary. All tracks required at 
least 14 hits in the 18 planes. The final fit at this point in the tracking routine gave 
the position of the track in x, y and z (called Xsms, Ysms and Zsms respectively) 
at the SM3 bend plane as well as the slopes TAN^a,(SM3) and TAN^j^(SM3), where 
TAN^-r and TAN^^^ are the ratios of the x and y momenta to the z momentum using 
the sign convention that TAN^-^ (or TAN^^^ ) is positive if or '^y is in the positive 
xory direction. 

The difference TAN^/(5M3) - TAN^^ ^{SM3), where U and D refer to up- 
stream or downstream of the SM3 bend plane, with the field map of SM3 determined 
the y-z projection of the momentum. Combined with TAN^j;(SM3), the four mo- 
mentum p*^) at Station 1 was fully determined. A final cut on the 18 plane fit 
requiring the track to have Xpdf ^^^^ than 5 was made. 

The full reconstruction of the track downstream of Station 1 was then com- 
pleted. Each track found above making all cuts was projected back to Station 4. 
Since there was considerable hadron absorber material throughout this projection, 
the projected intercept at Station 4 was compared to a 'window' of possible hits in 
the 3 proportional tube layers as well as the hodoscopes. These windows were wide 
enough to allow for multiple scattering within 5 standard deviations at each detec- 
tor plane around the projected straight line intercept. High momentum tracks had 
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windows of less than a few cells. The detectors were then scanned for hits within the 
windows. A track was further considered if it had 3 of the 5 planes (3 proportional 
tube planes and 2 hodoscope planes) registering a hit. 

This procedure was done for all candidate tracks found from the event on tape, 
and the resulting number of such tracks as well as numbers of hits, clusters and fitting 
information was passed to arrays for further use. 

3.3 Trace-back Through SM12 

Once all candidate tracks for an event were found, the event was passed to a section 
of code that would 'trace-back' each candidate track to the assumed place of origin on 
the z axis, Ztgt- For this Chapter, 'target' (such as tgt in Ztgt) is used to denote either 
the actual targets in the target wheel (target analysis), or the dump (dump analysis), 
unless it is necessary to distinguish between the two. This was done sequentially to 
each track in the order they were found in the tracking section of the code. Since 
there were no detectors upstream of Station 1, this was done using the magnetic 
field map of magnet SM12 as well as experimentally determined energy losses for the 
absorber materials [201; to 'project' the most probable path of the muon back toward 
its assumed origination. Once this projection was calculated, the effects of multiple 
scattering were added by use of a single bend plane approximation, resulting in the 
final calculated momenta for the muon at its assumed point of origin. 

The trace-back section of the analysis code is described here as it was originally 
written. The single muon analysis made several changes to the procedure which will 
be described in Section 13.5.11 however, the basic trace-back was preserved and valid 
for either dimuon or single muon events. Trace-back was done in three stages for each 
track found. 

The first stage entailed adding back in lost energy due to interactions within 
materials as well as changes in trajectory due to the magnetic field within magnet 
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SM12 beginning at Zsm3, the location of the bend plane in SMS. The trace-back 
procedure 'swam' the particle back to a fixed point along z, called Ztgt, which was 
either the location of the thin targets {z = —24.0 inches) for target events, or a 
fixed z location inside the dump for dump events {Ztgt for dump events was set to 
the distance at which 1/2 of the protons intering the dump would have interacted, 
Ztgt = 85.1 inches). 

The code began by adding in the beam offsets Xojf and Yojf, and then trace 
the muon upstream until it reached either the beginning of the copper dump at 2; = 68 
inches for target events, or until the next incremental distance in z (referred to as 
sections which were 2 inches in length) upstream of Ztgt where events from the dump 
were forced to originate from. 

Candidate tracks must pass cuts requiring the track to actually remain within 
the volume of the magnet during the trace-back. The x-y plane at Ztgt was called 
the analysis plane. The placement of the forced points of origination are discussed in 
detail in Section 13.4.21 

Corrections for energy loss used the mean energy loss for a muon in 2 inches of 
the material being traversed ISO]-^ The materials, in the order taken by the trace-back 
procedure, were 72 inches polyethylene, 108 inches carbon, 24 inches copper (called 
the copper wall) and either up to 168 inches of copper for target events or about 151 
inches of copper for dump events. The amount of the copper dump traversed was 
stored for use during the second stage of the trace-back. 

The second stage of the trace-back used the amount of dump material the 
track had traversed to determine the location along z at which correction for multiple 
scattering should be done. The analysis code used a single bend-plane approximation 
to correct for multiple scattering that occured while the muon traversed the dump and 
hadron absorber materials in magnet SM12. The single bend-plane approximation 

was based on determining a point along the z axis called Zgcat, where the effects of 

^The energy loss code calculated the mean energy loss for a muon of given momentum while 
travering the entire amount of absorber. The loss was then averaged over the 2 inch increments. 
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the multiple scattering could be approximated in a single large scatter, or bend (see 
for example |31j)- The location of Zscat was calculated using 

Zscat = 0' Idmp + 

where a is a percentage of the amount of dump material traversed and min is the 
minimum distance along z at which the plane could be set. Since all single muon 
events required full traversal of the dump, which restricted all hadrons to have a set 
open decay length, a was set to 0. 

The multiple scattering correction required retracing the muon from down- 
stream of magnet SM12 to Ztgt (at either the targets or in the dump), storing the 3 
momenta 'p^cat and position (X^^at, Yscat, Z^cat) found at Zscat during the trace-back, 
as well as the position {Xtgt, Ytgt, Ztgt) at Ztgt- ^tgt and Ytgt are commonly referred to 
as the uniterated x and y intercepts at the analysis plane. The difference between the 
uniterated intercepts and the beam centroids, AX and Ay , as well as an angular 
difference or correction to the track, AOr,. and A6y were calculated using 



(3.1) 





Ztgt — Zscat \ /\Y 



(3.2) 



The third stage was an iterative process of tracing the muon back to the anal- 
ysis plane from the scattering plane and testing the iterated x and y intercepts at 
the analysis plane against a distribution of acceptable values for Xtgt and Ytgt- If 
the current trace-back failed, new values for the angles AO^ and A6y were calculated 
and the process was repeated. The kinematical variables (all referred to as iterated 
variables) were calculated after this final step in the retracing of the muon. 
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3.4 Spectrometer Calibrations 



The analysis code required several input parameters to be set during initialization. 
The calibrations required were a scalar multiple for the magnetic field strength of 
SM12, called TWEE, the beam offsets Xofj and Ygfj, used to center the beam at Ztgt, 
the beam angles XSLP and YSLP and the position at which to place the scattering 
plane, Zscat, called ZSCPLN. 

Since the targets were relatively thin, it was assumed that the most accurate 
calibration for the magnetic field strength scalar, TWEE, would be found by fitting the 
dimuon data to the J/\E' mass (3.097 GeV c^^). The calibrations for dump analyses 
were then found using this field scaler, with a small variation (less than 2 percent) 
allowed for uncertainties in the field map as well as the effect of the thickness of the 
target and location of the plane of analysis used for dump events. All dimuon target 
analyses used data taken with the beam incident to all three targets while analyses 
for determining the dump calibrations only used data taken with the target wheel at 
the empty position. 

The calibrations were determined using an iterative process, changing the value 
for one of the variables and performing a new analysis of the Data Summary Tape. 
Results were plotted and subsequent changes were made to optimize the variable for 
that run. After all variables were determined for both target and dump the code was 
calibrated with the final settings and the initial reductions for single muons from that 
run were performed. Calibrations for this analysis are given in Table IB. 21 page 11471 

3.4.1 Cuts Used For Dimuon Events 

The cuts on dimuon events used for determining the optimal spectrometer calibrations 
were: 

1. Valid PhysAl, 2, 3 or 4 trigger must be set. 

2. Event must have only two candidate tracks. 
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3. Event must have two valid tracks after tracing section of code. 

4. Event must have oppositely charged muons. (Opposite sign pair.) 

5. The estimated z vertex, ZUNIN, must be within 50 inches of Ztgt. 

6. TAWy < 0.030 for both /i+ and 

7. TAm^ < 0.028 for both /i+ and /i". 

8. 2.1 < m^/, < 4.1 (GeV c^^). 

The cut on TAN^^ was to insure both muons from the event traversed the 
entire length of the dump as explained in Section 13.5.21 The cut on TAN6'2. was used 
to reduce the number of events where either muon may have scattered back into the 
acceptance after interacting in the walls of SM12, and m^^ is the reconstructed mass 
of the dimuon pair in GeV c~^. 

To select higher quality events, the E866 analysis code calculated an estimate 
of the position along the z axis where the parent hadron would have been created for 
the dimuon pair being analyzed. This calculation was based on information found 
as the two muons were traced back to the analysis plane at Ztgt during the second 
stage of the trace-back. The second trace-back section of code determined the four 
intercepts, X^t^ F^^^ as well as the two four momenta, (^, J)^) at Ztgt where ± 
refers to the two muons, /i+ and of the dimuon pair. The distance in y between 
the pair was calculated by F^^^ — 1^^^, and the estimated z vertex of the parent hadron, 
ZUNIN, was then calculated using: 



ZUNIN 



y- _ Y+ 

^tgt ^tgt 



(inches) (3.3) 



XTAWy- - TAN^/ 

The smaller the absolute value of ZUNIN calculated, the more likely the parent 
hadron was produced at or near Ztgt- 

The magnetic field strengths used in magnets SM12 and SM3 while taking the 
single muon data, plus the restriction that all events must have all muons trace back 
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Run 2753 Target Dimuons 
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Figure 3.2: Estimated z vertex, ZUNIN, for a target analysis of dimuon events from 
Run 2753 with any of the targets presented to the beam. Left is all events passing the 
dimuon cuts except for the ZUNIN and mass cuts, plotted as a function of ZUNIN. 
Right is same data with the mass cut, 2.1 < m^^ < 4.1 GeV c~^ , applied. Vertical 
lines indicate the limits on ZUNIN used in calibrating the spectrometer. Compare 
with figure EISl where no target was present. 

through the entire length of the dump, made the cut on the estimated z vertex of 
the dimuon pair, ZUNIN, insufficient by itself to isolate events originating from the 
target. Limiting contamination of dump events in a target analysis was accomplished 
using the cut on ZUNIN listed above plus the reconstructed mass cut. The effect of 
the mass cut on ZUNIN for the target analysis of run 2753 is shown in Figure 13.21 
Left is ZUNIN for all events with the beam incident on any of the three targets, and 
right is the same data with just the mass cut applied. Figure IX^ shows the reduction 
of dump events contaminating a target analysis using this cut by looking at data 
taken with the target position empty. Left is a plot of ZUNIN for a target analysis 
of run 2753 for all spills where the target position was empty. Right is the same data 
after the mass cut was applied. 



The absolute calibration of the luminosity monitor IC3 was unknown for these 
low intensities, but it was assumed that the total recorded by the monitor when the 
signal not busy was set (IC3SB) could be used to estimate the number of dimuon 
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Run 2753 Target Dimjons 
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Figure 3.3: Plots of the estimated z vertex, ZUNIN, for a target analysis of Run 2753 
where the target wheel was in the empty position. Left is the dimuon events passing 
all cuts except the mass and ZUNIN cuts. Right is the same data after the mass cut, 
2.1 < m^^ < 4.1 GeV c~^ , was applied. Vertical lines indicate the limits imposed on 
ZUNIN while determining the calibrations for the spectrometer. COmpare to figure 
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events in a target analysis that originated in the dump by comparing the appearant 
number of events in a target analysis when there was no target present. The number 
of events in a target analysis (which was the sum of the events from all the targets) 
can be estimated from: 

AT ^ 

where Nt{D) is the number of events in a target analysis that originated in the 
dump, A'o is the number of events in a target analysis performed for the empty target 
position, ICSSBi was the total IC3SB when target i was presented to the beam (the 
targets were labelled through 3 for empty, thin copper, beryllium and thick copper 
respectively) and Fi was the fraction of the beam incident on the dump: 

Fi = exp (^-k y-^ 

Here is the thickness of target i (cm), pi the density of the material in the target 
(gm cm~^) and Xj^i is the nuclear interaction length of the material of target i (gm 
cm~^). An example of the estimated dump contamination for Run 2751 is presented 
in Table EU (page . 

3.4.2 Calibrations 

The optimum settings for the beam offsets, Xofj and ^o//, as well as the magnetic 
field strength scalar TWEE were found using plots of the X^^ and Y^^^ intercepts at 
Ztgt for the /i^ individually, as well as the J/\l/. The analysis plane, Ztgt, was set to the 
known location of the targets, Ztgt = —24.0 inches. TWEE and Yoff were changed 
until the centroids of Yt'^t were centered about Ytgt = 0. Due to tracking limitations 
in the x direction, Xojf was varied until the centroids of X^t were equidistant from 
Xtgt = 0. Figure 1231 shows typical plots of X^t ^tgt the target dimuon analysis 
of run 2751. 
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Figure 3.4: Plots of Xtgt and Ytgt used in determining the beam offsets Xgff and Ygff. 
Figure is the final settings used for Run 2752. Vertical error bars are statistical only. 
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Beam angle corrections XSLP and YSLP (referred to as 6'^ and 6y below) 
were then determined for both target and dump analyses since the magnetic field in 
SM12 deflected the proton beam, using the two functions 



^Px _ 


'Px 










Pz 




Pz 





^Py _ 


'Py 




'Py 




^Pz 


Pz 


(/^+) + 


Pz 





The beam angle corrections were varied until the centroids of the two plotted 
variables were centered about 0. Figure shows these two variables plotted during 
the dimuon target analysis of run 2752. Figure 13.61 shows the reconstructed mass 
spectra for dimuon pairs calculated from target (left) and dump (right) analyses of 
run 2752 using the spectrometer settings used to create Figures 13.41 and 13.51 Table 
IB. 21 (page I147|) give the calibrations found for the runs used in this analysis. 



Run 2752 Target Dirruon 




-0.01 0.01 -0.01 0.01 



IP. /IP, IPy/lPz 

Figure 3.5: Plots of functions 13.41 and 13.51 used to calibrate the beam angles for Run 
2752. Vertical error bars are statistical only. 
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Run 2752 Dimuor Analyses 
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Figure 3.6: Reconstructed dimuon mass, m^^^ GeV c~^ for the target (left) and dump 
(right) of run 2752. Fits are to the J/^ (mj/^ = 3.097 GeV c~^). Errors are statistical 
only. Lines indicate the mass cut used in the dimuon analysis. 
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The optimal z position for the scattering plane, Z^cat-, was determined for 
target events by plotting the reconstructed dimuon mass, m^^, versus the estimated 
z vertex, ZUNIN. A sample of the method used for target dimuon events from Run 
2753 is shown in Figure 1X71 Top is m^^ versus ZUNIN with Zg^at = 150.0 inches, 
middle is the same data with Z^cat set to 175.0 inches and bottom is the result 
with Zscat = 200.0, inches which was found to give the optimum value for target 
events. Determination of the optimal value of Z^cat was found by iteratively changing 
the position, re-analyzing the dimuon data and fitting the reconstructed mass peak 
around the J/\E'. The optimal value was found by picking the position giving the 
minimum width. The optimal values were determined to be Zgcat = 200.0 inches 
for target events and Zg^at = 235.0 inches for dump events. Figure IXHl shows a plot 
of m^f^ versus ZUNIN for dump events from Run 2753 with Z^cat = 235.0 inches 
for comparison to the target analyses shown in Figure 13.71 The width of the mass 
spectrum for dump events was less sensitive to the location of Z^cat than target events 
due to the thickness of the proton interaction region and the limitation of having a 
fixed analysis plane. Figure 13.91 shows two plots of the reconstructed mass spectrum 
from run 2753 dump events. The left figure has Zgcat = 200.0 inches and the right 
Zscat = 235.0 inches. Both plots use the same calibrations otherwise. 



42 



Run 2753 Target DImuons 
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Figure 3.7: Reconstructed dimuon mass, m^^ (GeV c~^), versus the estimated z vertex 
(inches) for the target analysis of Run 2753. All plots use the same spectrometer 
calibrations except the position for the scattering bend plane, Zscat, which was set to 
150.0 inches for the top plot, 175.0 inches for the middle and 200.0 inches on bottom. 
Plots were used to determine the optimal position resulting in the minimum width of 
the mass spectrum, which occured when there was no slope. 
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Run ^753 Dump Dimuons 

ENTRIES 28401 



-60 -40 -20 20 40 60 

ZUNIN (inches) 

Figure 3.8: Reconstructed dimuon mass, m^^ (GeV c~^), versus the estimated z 
vertex (inches) for the dump analysis of Run 2753. Similar plots and fitted mass 
spectrum plots to those in Figure EiH were used to determine the optimal position of 
Zscat, resulting in the minimum width of the mass spectrum. 
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Figure 3.9: Reconstructed dimuon mass, m^^ (GeV c~^ ) for the dump analysis of 
Run 2753. Fits are to the J/"^ (w^j/* = 3.097 GeV c~^). Left is the reconstructed 
mass using Zgcat = 200.0 inches and right is the same data using Zgcat = 235.0 inches. 
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3.5 Single Muon Analysis 



Once the spectrometer had been cahbrated the analysis code was initiahzed for single 
muon reduction. This entailed setting several initial cuts to be made on all events as 
well as initializing different output to the n-tuple files. The single muon analysis was 
primarily interested in information regarding the tracking and trace-back, and effort 
was made to develop an estimate of the region along the z axis where the parent 
hadron originated. 

3.5.1 Single Muon Changes 

The original code used pairs of tracks that had been traced back. One criterion used 
to determine a single muon event was limiting the number of candidate tracks to 
one. The code was modified to allow single muon tracks to become valid events for 
the remaining sections of code by creating a ghost muon having the same ^ but 
reversing the direction of "p^y. All single muon events thus became opposite sign 
dimuon events for analysis in the remaining sections of code, with the exception that 
ZUNIN = 0. Since all dimuon information was output for a muon pair (usually 
opposite in charge), the single muon code was changed to differentiate which track 
was the real track and what charge that muon had. This was required since the y 
dependent ghost track information was reversed in sign and placed into the output 
n-tuple as the opposite sign muon of that pair, and the original code was developed 
for only one current direction in SM12, creating an ambiguity in the analyzed events 
as to which muon was really the fi^ and which was the fi^. Since the original code 
reconstructed events using both muons, this ambiguity had no effect on their results 

A short section of code, after the pair of muons had been fully traced back, 
allowed for rotation of the z axis to allow for beam angle corrections. The new 
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momenta were calculated from: 

bLI = cos(C) \-pA+sin{d'^) 1^,1 (3.6) 

\p'y\=cos{e'y) \-fy\ + sin{e'y) 1^,1 (3.7) 

b'J = ^J\'p.V-\p'^^-\p'y\^ (3.8) 

where primed implies after rotation. The angle corrections, 9'^ and 9'y were found using 
a large number of dimuon events for each run as described in the previous section. 
Figure EUni shows the ratio {p^ —pt)/pt versus pt for /i+ (top) and /i~ (bottom) single 
muon events from a target analysis of run 2755. 

The original code was also altered to allow calculating the 'z of closest ap- 
proach' for three cases. The code was changed to calculate the distance a;, y and 
R = \J + y"^ from the muon to the z axis at each section during the first trace- 
back. The code then output the z locations at which each of these became a minimum. 
These minima gave a qualitative examination of the ability of the tracing section of 
the code to determine the origination of the muons in z, i.e. target or dump, for 
single muon events. 
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Run 2755 Target Single Muons 
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Figure 3.10: Effects of small angle rotation correction applied to single muon events 
from the target analysis of Run 2755. Events have the trigger, tracking, circle and 
maximum tan{9y) cuts applied. Top is /i'^ and bottom is ii~ . 
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3.5.2 Single Muon Cuts 

Initial analysis of the data for single muon events used a series of cuts designed to 
insure that the event was the result of a leptonic or semi-leptonic decay resulting in 
one muon, the decay muon must have traversed the full length of the solid copper 
dump and the event must have originated from the correct region along z. Triggers 
and number of candidate tracks for the tracking section of the code were used for 
the first part, and cuts developed using both iterated as well as uniterated physics or 
tracking variables were used for the others. Subsequent analysis used cuts designed 
to further restrict these requirements as well as ensure the validity of the tracking 
section of the analysis code. 

Trigger Cut 

A cut was placed after the trigger was read out and any event not having a PhysBl 
trigger set was cut. A further restriction on which triggers were not allowed was 
placed so that any event having any of the five dimuon triggers consisting of PhysAl, 
PhysA2, PhysAS, PhysA4 and Diag3 was cut. This combination is referred to as 
the trigger cut. This cut limited all events to those having the required X134L/R 
coincidence and having no second track setting a dimuon trigger. 

One Track Cut 

Once an event had been read in and passed the trigger cut, the code began construct- 
ing all candidate tracks from hit clusters in Stations 1, 2 and 3 for that event. Once 
all candidate tracks had been found for that event the code cut any event having more 
than one candidate track. A second cut was placed after the event was passed from 
the retracing section of code to insure no event had more than one retraced muon. 



48 



Fixed Tracking and Retracing Cuts 

Any muon that tracked outside the apertures of any wire plane, hodoscope layer or 
the volume in magnet SMS was cut during the tracking section of code. During the 
trace-back section of code the analysis cut all events having x and y slopes, TAN6'a. and 
TAN6'j^, pointing outside the physical volume inside magnet SM12. These physical 
apertures were checked continuously during the trace-back section of code for each 
muon. 

Projection Cuts 

To insure that the single track could have actually set the X134L/R coincidence 
trigger from the hit clusters used in the tracking section of code, the Data Summary 
Tapes were subsequently analyzed with a cut imposed on the minimum projected 
distance in x, referred to as Xsi, from the z axis in the x — z plane at z = 770.72 
inches and the maximum distance in x, referred to as Xs3, &X z = 1822.0 inches. The 
two z positions are the two hodoscope banks XI and X3 respectively (see figure E31 
1221 for a visual explanation). Hodoscope plane X\ set the minimum opening angle in 
X that a muon must have to set the X134L/R coincidence, and hodoscope plane X3 
set the maximum opening angle the muon could have to set the same coincidence. 
This cut not only insured that the track would set the trigger, but when combined 
with the circle cut described in a later section, they also cut events that had scattered 
back into the acceptance from the walls of the spectrometer magnet. 

Figure IS'. 1 II shows the effect on TAN^^. when these cuts were applied to /i^ from 
the initial dump analysis of Run 2748. An event satisfying the X134L/R trigger cut 
should, with no multiple scattering, have a minimum |TAN^a;| greater than 0.017, and 
a maximum of 0.028. Data outside of these limits are events that, while traversing the 
dump and absorber materials, either scattered into the maximum acceptance angle 
(|TAN^a;| > 0.028) before hitting the spectrometer or scattered into the minimum 
angle (|TAN6'2,' | < 0.017). Similar effects were seen from the target. 
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Circle Cut 

A cut on the maximum allowable \Xtgt\ and \Ytgt\ was used to begin the separation 
of target and dump events in the initial target and dump reductions. This initial 
cut was imposed after inspection showed that events originating from Zfgt primarily 
retraced (in a target analysis) to a relatively small distance from the z axis when 
plotted as Ytgt versus Xtgt at Z = Ztgt- Figure EHH illustrates this effect on the data 
when the analysis plane is set to the target location. The left plot is events when any 
of the three targets were incident to the beam during run 2748. The plot is limited 
to events where both Xtgt and Ytgt are less than 10 inches away from the z axis. The 
data has the trigger and one track cuts applied. In contrast, the plot on the right 
is the same data where the target position was empty, which means that all of the 
events originated from hadrons created in the dump. 

The initial reduction cut all events having the magnitudes of either Ytgt or 
Xtgt greater than 2 inches. The final reduction used a cut on the maximum distance 



TAN6'j,and Momentum Cuts 

To limit the open decay distance to one value the single muon analysis restricted 
all muons to traverse the entire length of the copper beam dump. There were two 
original methods under consideration for this cut. The first was to cut any events 
that exited the dump as they were traced back to Ztgt- However, the acceptance for 
dump events would have been considerably larger since, for the most part, all events 
originating from hadrons produced in the dump traversed the entire length left over. 
That was not the case for target events. 



defined as D 




Xtt + Ytqf Events having D > 1.0 inches were cut. 
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Figure 3.11: Effects on the slope tan{6j) with cuts apphed on the minimum distance 
for Xsi (upper right), maximum distance for (lower left) and both with the circle 
cut added (lower right). Upper left is all //"^ for single muon events from the initial 
target analysis of Run 2748. 
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Run 2748 Single Muon Target Analysis 
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Figure 3.12: Ytgt vs Xtgt for data from Run 2748 analyzed at Ztgt — —24.0 inches. 
Left plot shows data taken with any target, right is the data with the target position 
empty from the same run. 
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The cut used consisted of determining the maximum and minimum slopes in 
the y — z plane combined with a minimum momentum muons must have and still 
traverse the length of the dump. The same cuts would then be applied for muons 
originating from either the targets or the dump. Studies of target events using the 
traced y position of the muon at z = 86 inches (referred to as y86) were performed 
using a variety of cuts in both TAN^^and momentum, I'p'l. It was decided the cut 
should remove the fewest high transverse momentum events possible and still meet 
the criteria of forcing all muons from hadrons produced in the target to retrace the full 
length of the dump. The cuts selected for data taken with the spectrometer magnets 
having parallel fields were |TAN^j/| < 0.030 and \~p \ > 55.0 GeV/c. 

Figure IH. 1 31 shows the effects of these cuts. The top two plots show the from 
a target analysis of Run 2748 having the trigger, tracking and circle cuts applied on 
the left, and the additional TAN^^and I'p'l cuts applied on the right. The bottom 
plots are the same data for the 

The two figures on the left show the effects of the increased decay length given 
light hadrons that may pass under or over the copper beam dump. The decay length 
for events that must traverse the entire length of the dump is 92 inches, while for 
those passing above or below it is 259 inches. The decay muons from these hadrons 
are those in the smaller peaks having larger mean distances from Y8Q = 0. 

Dimuon data had suggested that magnet SM12 had dropped approximately 
0.40 inches in the y direction. The bottom of the dump would therefore heaty = —3.9 
inches, and the top at y = 3.1 inches. The application of the TAN6'j, and cuts 
together force all decay muons from hadrons originating from the target to trace 
between those values at Z = 86 inches. 
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Run 2748 Single Muon Target Analysis 
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Figure 3.13: The values of y86 for //+ (top figures) and fj,~ (bottom figures) from a 

target analysis of target events in Run 2748. Left are events with the trigger, one 
track, aperture, projection and circle cuts applied, right are the same events with the 
TAN9y and minimum | P \ cuts added. The two figures on the left show the effects 
of the increased decay length given light hadrons that may pass under or over the 
copper beam dump. The decay length for events that must traverse the entire length 
of the dump is 92 inches, while for those passing above or below it is 259 inches. The 
decay muons from these hadrons are those in the smaller peaks having larger mean 
distances from Y86 = 0. 
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3.5.3 Single Muon Cuts Summary 

The single muon analysis has initial cuts as well as some intermediary and final cuts. 
These are: 

• Initial Cuts: 

1. PhysAl trigger requirement. 

2. No dimuon trigger set. 

3. One candidate track. 

4. One traceable track. 

5. Spectrometer aperture cuts, generally; 

1^1 > 30.0 GeV/c. 
|TAN^^ I < 0.040. 
\TAN9y I < 0.050. 

6. Loose radius or circle cut set to: 



• Intermediate Cuts: 

1. Loose projection cuts set to: 

Xsi > 11.0 inches. 
Xs3 < 26.0 inches. 

• Final Cuts: 

1. Tight radius or circle cut set to: 



D^JXL + Y,l,< 2.0 inches. 





2. Tight projection cuts set to: 
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Xsi > 12.0 inches. 
Xs3 < 25.5 inches. 

3. Minimum momentum set to: 

\^\ > 55.0 GeV/c. 

4. Maximum y slope set to: 

\TAN9y I < 0.030 

3.6 Data Reductions 

The analysis required four separate initial reductions for each tape, one each for 
target and dump dimuon and one each for target and dump single muon events. The 
results from each reduction were stored as Data Summary Tapes, or DSTs. The 
basic procedure was to first reduce the data for dimuon events from the targets, 
then dimuon events from the dump. These DSTs were then used to calibrate the 
spectrometer, after which the data tapes were reduced to the two initial single muon 
DSTs. 

3.6.1 Final Analysis and Presentation of Results 

The final analysis of the data was performed using the CERNLIB software Physics 
Analysis Workstation, or PAW. This software is available from CERN for most plat- 
forms |17j. The final output from the analysis code was in four forms, one a de- 
tailed log file of the analyses performed and relevant information to the run, a Data 
Summary Tape for successive analyses, an analysis generated set of histograms for 
variables of interest under several conditional criteria such as cuts and the charge of 
the muon using the CERNLIB routine HBOOK, and an n-tuple file containing all 
events that are of interest with the desired information for each of those events as a 
database. 



56 



The information in each n-tuple file was searched for those events of interest 
and placed into histograms. The histograms could be defined for any binning width 
within reason, but the final histograms for the transverse momentum had 40 bins 
from to 10 GeV/c resulting in a bin width of 0.25 GeV/c. Reconstructed masses 
from dimuon events typically had bin widths of 0.60 GeV c~^. The use of PAW and 
n-tuple files greatly increased the ability to change search criteria without lengthy 
re-running of the analysis code. 

The single muon analysis increased the original number of variables in the n- 
tuple file for each event to accommodate increased tracking information as well as 
uniterated momenta and momenta before rotation by the beam angle corrections. 
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Chapter 4 
Monte Carlo 



Central to determining the open charm cross sections was the use of the E866 Monte 
Carlo to produce expected muon decay spectra from the hadrons that contribute. 
The hadrons that contribute significantly to the single muon spectra are vr 



to the production of and / D since the spectrometer could not differentiate 
between these open charm mesons and any other open charm or heavier mesons such 
as the D* or Ds since those decay strongly to one of the four open charm mesons 
used. 

Muons resulting from the decay of mesons with large momentum usually have 
smaller transverse momentum than the parent. This is due to the fact that the 
momentum of the parent is 'shared' between the decay products, and pr is a function 
of the momentum of the particle, px = |'p'|sin(^). Figure Wl] shows the resulting 
Pt (/i) distribution from open charm hadrons thrown with p^ (h) between 8.00 and 
8.25 GeV/c . 

For brevity, and to avoid confusion, the transverse momentum shift will be 
referred to as the pt shift, and the transverse momentum of hadrons will be denoted 
Pt (h) while the transverse momentum of muon will be denoted pt (/i). Also, since 
there are two charges of muons, as well as hadron/anti-hadron, the usual reference 




The open charm cross section was limited 
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will be to hadrons and muons unless the distinction between hadron/anti-hadron or 
/i+ and /i~ is required. 

The pt shift imposed the requirement that all hadrons needed to be thrown to 
larger pt {h) than the maximum pt (/i) in the data. Hadron production falls steeply 
as a function of pt {h) (see Equation 14.21 fpage and Figures and lOl on pages 
EniandEZI) The value of the differential cross section, which falls exponentially or 
nearly so as a function of pt-, is a measure of the number of hadrons that will be 
produced. Worse, the likelihood that a hadron would decay also falls exponentially 
as a function of the momentum of the hadron. The minimum momentum of the 
muons from the data after final analysis was 55 GeV/c, and since the momentum 
is shared between the decay products, the minimum momentum for hadrons would 
necessarily have to be greater than the minimum for the muons. The decision was 
made to throw the hadrons flat in transverse momentum, center of mass rapidity y 
[y = \ log j and 0, and then use weights to shape the resulting spectra. 

Use of the term fiat means that a random value was thrown between a lower and 
upper limit, so the distribution was 'flat' when it was histogrammed. 

The E866 Monte Carlo had two distinct parts in its original form, an event 
generator section and a decay muon tracing section for dimuon events. Usual prac- 
tice was to use an outside source of dimuon events generated in other custom event 
generators using hadron production based on programs such as PYTHIA, where the 
resulting dimuon kinematics were read directly into the tracing section of the Monte 
Carlo. The muon tracing section of the Monte Carlo was used to simulate the E866 
spectrometer and produce simulated tracking information from the four tracking sta- 
tions as well as trigger information from the hodoscopes. The output trigger and 
track information was placed into an E866 data format file for later analysis. Other 
useful information, such as the weights used in this analysis, could be passed to the 
analysis routine and read out for each event. The information could then be used 
directly in the analysis and/or written to the n-tuples if desired. 
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Several major modifications were made to the Monte Carlo to simulate single 
muon decays from the hadrons of interest. The event generator was modified to throw 
the desired hadrons using a distribution fiat in px (h), center of mass rapidity, y, and 
(p. Once the hadron was thrown the event generator was required to determine a decay 
point based on the decay properties of the hadron and its momentum. Traversal of 
target or dump materials and magnetic field effects of the spectrometer magnet SM12 
were simulated in a new section of the program for all charged hadrons. Hadron decay 
was performed using a CERN routine named GENBOD, and all decay products were 
transformed to the lab frame using the CERN routine LORENB. Secondary hadron 
decays, when required, were introduced by repeating the hadron tracing and decay 
sections for each secondary hadron of interest resulting from the primary hadron 
decay. These modifications allowed for large numbers of events at high transverse 
momentum since all hadrons that made it into the spectrometer were allowed to 
decay, and all hadrons were thrown fiat in pt (h). Figure shows a representation 
of the single muon Monte Carlo process. 

Various cuts on primary hadrons, secondary hadrons and decay muons were 
implemented throughout the event generator and muon tracing sections of the Monte 
Carlo. The first set of cuts were loose aperture and momentum cuts taken on the 
hadrons immediately after they were thrown and histogrammed. This increased the 
efficiency of the Monte Carlo since similar or harder cuts on the thrown hadron were 
done later while the hadron was being traced through the spectrometer magnet. All 
decay muons or secondary hadrons having \ 'p'\ < 50 GeV/ c were cut immediately after 
the decay and boost routines were done. All hadrons that were allowed to contribute 
secondary muons that resulted in more than one muon traversing the spectrometer 
were cut. A cut was imposed on all muons that tracked through any of the pulled x 
hodoscopes. 

After a muon had passed all cuts the parent hadron and decay muon informa- 
tion required for producing weighted spectra was written to the MC data file. 
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Figure 4.1: The muon distribution from open charm hadrons thrown with random 
transverse momentum between 8.00 and 8.25 GeV/c. 
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Figure 4.2: Basic outline of the Monte Carlo process to calculate the various contri- 
butions to the single muon spectra. 
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4.1 Weighted Spectra 



Extraction of the open charm cross sections required a set of Pt (/i) spectra from all 
hadrons contributing to the data. These spectra were required to be the muon spectra 
that would result for a given number of protons interacting in the various materials 
of the targets and copper dump. A set of weights were given each muon as it was 
written to disk. These weights are shown in Table 14.11 The total weight of a muon is 
the product of the individual weights (listed in Table l^!T|) appropriate for the event. 

Table 4.1: Weights calculated for muons in the Monte Carlo spectra. Hadrons are 
divided into groups, light (L), open charm (H) and both (B). {mbc^ GeV-^) 

is the p-p differential cross section at the thrown p^ (h) and cm. y of the hadron, 
j[a{pT(h)) jg ^Yie nuclear dependency of the hadron, Pd is the probability that the 
hadron has decayed in the spectrometer, P/v^ is the likelihood that a hadron produced 
in the target did not interact in any remaining material in the target and br{mode) 
is the muonic branching fraction of the decay mode. 



Weight 


Hadron 


Wa 
Wd 
Wna 


E^ {mbc^GeV-^) 

Pd 

Pna 
br {mode) 


L 
L 
B 
B 
B 



4.1.1 The Light Hadron p-p Cross Section Weight 

All hadron distributions were parameterized using p-p differential cross sections, and 
the following experimentally determined relationship was used to convert those cross 
sections to p-A cross sections: 

E-i^ = ^"('■"•'X E (mb GeV-) (4.1) 

dp-^ dp-^ 
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where h is the hadron of interest and ^^"(^'PtCi)) jg referred to as the nuclear depen- 
dency term. 

A parameterization of the p-p differential vr^ and K"^ cross sections as a 
function of transverse momentum, pT{h), and cm. rapidity y, developed by the 
British-Scandinavian Collaboration (BSC) at CERN ^| was used to calculate the 
cross section weight, (mb c'^ GeV~^) . The parameterization of the differential 
cross sections developed by the BSC is 

^ _ ptWY 

E '^^^{pr (h) ,y)=A, e^--^'^)) e(--^^) + (4.2) 
«r [pt {hy + M2) 

where values for the parameters Ai, B, D, A2, M, m and n suggested for use by 

the BSC, their best central value and undertainty are listed in Table 14.21 It should 

be noted that the exponent n was held fixed at 4 for all fits, and the collaboration 

suggested no deviation for that parameter. The center-of-mass (cm.) momentum of 

the proton, p^^^, introduces the energy dependency of hadron production and for 

this experiment was set to p'^am — 19.37 GeV/c. 

It was assumed that <y{K^) + cr(K^) = a{K^) + a{K^) and that a similar 
production relationship existed between the K'^ and as between the and 
K~ due to conservation of isospin. This was introduced into the Monte Carlo by using 
(t(X°) = a{K+) and a(K^) = a{K-). The Monte Carlo did not throw /K^ but 
instead threw Kl where the parent was either a K'^ (50%) or 1^ (50%). 

Figures (4.31 and (4.41 show the vr^ and cross sections for the parameteriza- 
tions in at ?/ = (cm.) for the energy of this experiment, ^/s = 38.8 GeV. 

4.1.2 Nuclear Dependency Term 

It is well known that the proton-nucleus (p-A) cross section, at least for light hadrons, 
does not scale to the proton-proton {p-p ) cross section by the number of nucleons, 
A, in the target nucleus. The Chicago-Princeton Collaboration (C-P) studied the 
production of light hadrons on various nuclear targets and developed what has become 
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Table 4.2: Values of the parameters used in the BSC parameterization of the vr^ and 
differential cross sections. Fits were to data with pt between and approximately 
6 GeV/c. Taken from |32j . 





Ai 


B 


D 


A2 


M 


m 


7r+ 


210 ±4 


7.58 ±0.11 


0.20 ±0.01 


10.7 ±0.7 


1.03 ±0.03 


10.9 ±0.4 




205 ±4 


7.44 ± 0.09 


0.21 ±0.01 


12.8 ±0.9 


1.08 ±0.02 


13.1 ±0.4 


K+ 


14.3 ±0.4 


6.78 ±0.21 


1.5 ±0.1 


8.0 ± 1.1 


1.29 ±0.03 


12.1 ±0.8 


K- 


13.4 ± 1.0 


6.51 ±0.23 


1.8 ±0.1 


9.8 ± 1.7 


1.39 ±0.04 


17.4 ± 1.0 
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Pt (GeV/c) 

Figure 4.3: The n'^ (red) and vr^ (blue) differential cross section as used in the Monte 
Carlo. Distribution is tak;en at cm. y = and = 19.37 GeV/c. The value of 

the differential cross section at the thrown px {h) and y was calculated and passed to 
the output MC data file and later used as a weight (Wx) for the event. 
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Figure 4.4: The (red) and K~ (blue) differential cross sections as used in the 
Monte Carlo. Distribution is taken at cm. y = and p'^am — 19.37 GeV/c. The 
value of the differential cross section for the thrown pj- [h) and y was calculated for 
each event and passed to the output MC data file and later used as a weight iWx) 
for the event. 
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known as the nuclear dependency term, relating the p-A to p-p cross sections by the 
relation 

where A is the nuclear number (referred to as the atomic weight in [SS! and h is the 
species of hadron. 

Light hadrons thrown in the Monte Carlo had the nuclear dependency cal- 
culated as a weight, Wa, that was passed to the Monte Carlo data files for use in 
weighting the thrown spectra after analysis. Since the open charm cross sections 
were determined as the p-Cu and p-Be cross sections, no nuclear dependency for 
the production of open charm was used. 

Data points from the Chicago-Princeton Collaboration (C-P) study for the 
nuclear number (or equivalently the atomic- weight) dependency for charged pion and 
kaon production from 400 GeV/c p-A collisions [34j were fitted to simple quadratic 
equations of the form pi + P2PT + PsPt using the CERN program MINUIT. The 
resulting parameterizations were used to calculate the value of the dependency for 
each light hadron thrown. Figure 14.51 shows the four fits to the C-P data, and Table 
14.31 gives the parameter values of the fits. The K~ parameterization was fitted to two 
regions in pj- {h). 

4.1.3 Probability To Decay 

All hadrons thrown were given a weight, Wd, based on the likelihood that the hadron 
would decay, L^), in the spectrometer. Since all hadrons must traverse some material, 
the loss of hadrons due to catastrophic collisions was included in the decay likelihood. 

The probability for a hadron with momentum | to decay in a distance z is 
given by 

Pd = 1 - exp (^-z (4.3) 

where z is the distance the hadron has travelled (cm), m is the mass of the hadron 
(GeV/c^), p is the momentum (GeV/c) and r is the proper lifetime (sec). The 
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Table 4.3: Values for parameters of the equation a{h,pT (h)) = pi + P2PT + PsPt- 
The K'^ used the same parameter values as the and the K used the same values 
as the K~. The K~ fit required using two regions; hadrons thrown with transverse 
momentum less than 4.0 GeV/c, and hadrons having transverse momentum greater 
than 4.0 GeV/c. 





Pi 


P2 


P3 


7r+ 


0.8074 ± 0.0072 


0.1368 ± 0.0060 


-0.0148 ± 0.0010 


7r~ 


0.8043 ± 0.0059 


0.1412 ± 0.0040 


-0.0146 ± 0.0005 


K+ 


0.8464 ± 0.0154 


0.1380 ± 0.0113 


-0.01481 ± 0.0019 


K- 


0.9295 ± 0.0148 


0.0477 ± 0.0121 


0.0053 ± 0.0006 




0.3649 ± 0.0654 


0.3289 ± 0.0266 


-0.0300 ± 0.0032 
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probability for a hadron to interact after traversing a distance z in a material is 
calculated from 

Pi = l- exp (^-^Y^ (4-4) 

where p is the density of the material being traversed (gm cm~^) and A/ (gm cm~^) 
is the nuclear interaction length of the material for the hadron, which is slightly 
dependent on the hadron momentum. The two processes are independent of each 
other, so the likelihood that a hadron has neither interacted nor decayed in a length 
z and will decay in the next differential thickness dz is given by 



dz 



m 
— exp 
pr 



m p 
pr \i 



(4.5) 



The fraction of hadrons that will decay while traversing a material is found 
by integrating Equation 14.51 and normalizing the result such that the likelihood for 
a hadron to either decay or interact while traversing a target of infinite thickness is 
unity. The fraction of hadrons that will decay while traversing a material of thickness 
z is then given by 



m A/ 



exp 



m Xj + ppr 

while the fraction that will interact, Fj, is given by 



m p 

— + f 
pr Xi 



(4.6) 



Fi 



PPT 



exp 



—z 



m p 

— + f 
pr Xi 



(4.7) 

l — Fj was 



mXi + ppT 

The fraction of hadrons not interacting while traversing a material, F^i 
more useful in calculating the likelihood for decay in the spectrometer. 

The assumption was made that all hadrons traversing the dump would either 
interact or decay, so the length of the dump was set to infinity in Equation 14.61 for 
simplicity. This assumption reduced Equation 14.61 to 

_ mXi 

i'D - 7 \ 

m Xi + ppT 

. It was also assumed that all open charm hadrons thrown from the target would 
decay if they didn't interact while traversing any remaining material in the target. 
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The likelihood that the hadron thrown would decay in the spectrometer was 
then determined as follows: 

• If a light hadron was thrown in the target, the likelihood for decay was calculated 
using 

Wd = F^jPn + (1 - FniPd) Fd 

• If an open charm hadron was thrown in the target, then 

Wd = Fmi 

• All hadrons thrown in the dump used 

Wd = Fd 

In all cases, the value of F/v/ used the difference between the thickness of the 
target and the z location where the hadron was produced f Section 14. 1.41 below), and 
the distance used in calculating Pr, was set to 92 inches (233.68 cm), the distance 
between the targets and the beam dump. 

The nuclear interaction length. A/, for light hadrons was calculated using the 
results published by A. S. Carrol, et al. [33]. The code initialized three interaction 
lengths associated with three momentum ranges of the hadron, < < 130, 
130 < Pt ^ 240 and pt > 240 GeV/c. The nuclear interaction length, calculated 
using the absorption cross sections determined by A. S. Carroll, et al, are given 
in Table 14.41 The nuclear interaction lengths for open charm was calculated by 
decreasing the absorption cross section for strange particles by 10 percent. 

4.1.4 Proton Interaction And Hadron Decay Distributions 

Two routines were used to model the distribution of hadron interactions in the mate- 
rials; one using a two dimensional Gaussian distribution to model the axial dispersion 
of the proton beam, and one using a distribution based on the likelihood that a hadron 
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Table 4.4: Absorption cross section, a a-, and nuclear interaction length, A/, for hadrons 
thrown in the Monte Carlo. All hadron momenta are in GeV/c. Nuclear interaction 

length for neutral kaons and / D are the same as their charged counterparts. The 
density for beryllium used was 1.848 gm cm~^, and for copper 8.96 gm cm~^ . The 
nuclear interaction lengths for protons was taken directly from the tables provided 
by the Particle Data Group. 



H pt 

Range 


Copper 
(Ta (mb) A/ gm cm~^ 


Beryllium 
(7a (mb) A/ gm cm"^ 


7r+ - 130 
130 - 240 
> 240 


611.6 ±34.8 172.5 ±9.8 
616.6 ±35.2 171.0 ±9.8 
608.6 ± 37.0 173.4 ± 10.5 


137.5 ±6.2 108.8 ±4.9 

139.3 ±6.2 107.4 ±4.8 

140.4 ±6.6 106.6 ±5.0 


Tr~ 


611.8 ±34.7 172.5 ±9.8 
611.4 ±34.7 172.6 ±9.8 
610.6 ±36.6 172.8 ±10.4 


138.1 ±6.1 108.4 ±4.8 
138.8 ±6.2 107.8 ±4.8 
139.7 ±6.4 107.1 ±4.9 


K+ 


543.5 ±31.1 194.1 ±11.1 

562.7 ±32.3 187.5 ± 10.8 

555.8 ±36.8 189.9 ± 12.6 


115.9 ±5.2 129.1 ±5.8 

122.2 ±5.5 122.5 ±5.6 

122.3 ±6.2 122.3 ±6.2 


K- 


572.7 ±33.0 184.2 ±10.6 
571.3 ±32.8 184.7 ±10.6 
570.0 ±35.6 185.1 ±11.6 


125.1 ±5.7 119.7 ±5.5 

125.2 ±5.7 119.5 ±5.4 
126.4 ±6.2 118.4 ±5.8 


D+ 


215.7 ±12.3 
208.4 ± 12.0 
211.0 ± 14.0 


143.4 ±6.5 
136.1 ±6.2 
135.9 ±6.9 


D~ 


204.7 ± 11.8 
205.2 ±11.8 
205.7 ±12.9 


133.0 ±6.1 
132.8 ±6.0 
131.5 ±6.4 


P 


134.9 


75.2 
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had not interacted in a randomly chosen distance. As before, the distance along the z 
axis was used as the total distance ignoring any small x oi y components. The beam 
distribution was provided in the original Monte Carlo. All lengths and distances 
throughout this section are in cm. 

All hadrons thrown assumed that a proton-nucleon interaction had occured, 
and that all hadrons would decay (the decay weight, Wo, described in 14.1.31 above 
provided the proper scaling for loss due to interactions). The Monte Carlo used an 
iterative process where a trial distance in the material at which the interaction or 
decay would occur (A^, and then tested the trial distance against the interaction 
or decay distribution. If the trial distance fell within the distribution, the routine 
then determined the point (A^;, A^,, A^) at which the decay or interaction would have 
occured based on either the distributed beam diffusion and beam correction angles, 
or the momentum of the hadron "j?. If the trial distance failed to fall within the 
distribution, the routine initiated a new trial A^ until a distance was accepted. 

Specifically, to determine the point along the z axis at which the proton had 
interacted for hadrons thrown from the target, the routine threw a random number (all 
random numbers were between and 1) and multiplied the number to the thickness 
of the target, providing the trial interaction distance, A^. The likelihood that the 
proton had survived to the trial distance was calculated using: 



where pA (gm cm~^) is the density of the target material A and A/(p, A) (gm cm~^) is 
the nuclear interaction length for 800 GeV/c protons in the material. L^j was referred 
to as the test function. A second random number, rn was thrown and compared to 
the value Liv/. If rn < L^i the trial distance A^ was accepted and the routine then 
determined the interaction point (X/,1/,Z'/) where Xi = 6^ + X SLP A^. 6^ was 
the distance in x away from the z axis determined by the beam dispersion routine, 
and XSLP is the beam angle in the x — z plane defined in 13.4.21 If rn > L^v/, the 
process was repeated until the test distance A^, was accepted. Yj was calculated by 
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substituting y for x, and Zi — —60.96 + A^. 

For hadrons thrown in the dump, the process to determine where the proton 

interacted was the same, except the trial distance was hmited to the distance at which 
95 percent of the proton interactions would have occured, and magnetic field effects 
were included. The maximum distance A^"^ was calculated using 



The trial distance was then calculated using A^ = rn A^"^. Once a z decay distance 
had been found, the routine projected the proton through the spectrometer to the 
point Zi = 203.2 + A^. The proton was given initial 6^ and 6y components as well as 
the correct beam angle before the projection was performed, so the interaction point, 
{Xi, Yi, Zi), was the values of the position of the proton at Zj after the projection. 

The distribution of hadron decays for hadrons produced in the dump were 
handled in the much the same manner. The test function used was 



where m is the mass (GeV c^^), p the magnitude of the momentum "p^(GeV/c), r 
the proper lifetime (s) and Xi{h, Cu) the nuclear interaction length (gm cm~^) of 
the hadron h. The maximum allowable distance, A^"^ was determined by setting 
Lni — 0.05 and solving for z. 

Two test functions were required for light hadrons thrown from the dump, one 
for the distribution of decays in the open decay region between the target and the 
dump, and one for the distribution of decays in the dump if the hadron survived the 
open decay length. The test function used for A^ < 233.68 (cm) was: 



where odl signifies the open decay length (233.68 cm). The test function for 233.68 < 

A^ < A™""^ (cm) was given as 



max 



/n(0.05) 



\i{ip, Cu) 



PCu 
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where P^"^' = 1 — exp [(—233.68 m) / (pr)] (the probabihty the hadron did not decay 
in the open decay region) and A™"^ was calculated by setting L^i = 0.05 in the 
second test function and solving for z and adding 233.68 (cm). Figure shows a 
plot of the two test functions used to determine the decay distribution of charged 
kaons having I'p'l = 55.0 GeV/c. 

1-- 

0.8: 

E 0.6- 
o 

5 0.4: 

0.2- 



^ 50 100 150 200 250 300 350 

cm 

Figure 4.6: The two test functions used to determine the decay distribution for 55.0 
GeV/c charged kaons thrown from a target. The value of the test function at the 
randomly chosen decay distance was compared to a random number between and 
1. If the random number was smaller than or equal to the value of the test function, 
the decay distance was accepted, if not a new decay distance was calculated and 
the process repeated until a distance met the required test. Distance is the distance 
downstream of the target {Ztgt = —60.96 cm) 

4.1.5 Decay Modes And Branching Fractions 

The Monte Carlo used a set of decay mode routines to determine the mode and 
branching fraction of the mode chosen. Since the Monte Carlo spectra were weighted 
spectra, only semi-muonic decay modes were allowed for open charm hadrons. Table 
l4.5l shows the decay modes allowed for each hadron species thrown in the Monte Carlo. 
The weight given each decay was the branching fraction listed. Strange hadrons were 
allowed to decay to charged pions, which were subsequently passed through the decay 
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routine as secondary decays. However, if more than one muon was accepted the event 
was cut due to the requirement in the analysis section where only one valid track was 
allowed for each event. 
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Table 4.5: The decay modes allowed in the Monte Carlo. The weight Wj^r given any 
muon that was accepted was the branching fraction as shown. Only the hadron decay 
modes are shown, the anti-hadron modes are the charge conjugate of the mode listed. 
The Kl hadrons were assumed to be from 50 percent of the time, and K 50 
percent of the time. Decay modes of the anti-hadrons are the charge conjugate of 
those shown. 
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Chapter 5 

Method Used To Extract The 
Open Charm Cross Sections 

This analysis used the single muon data and Monte Carlo muon spectra to extract 
four inclusive open charm differential cross sections, E — '^^g^ ' {pr), E — j^—^ (pt), 
E '^^ "^J^'' (pt) and E ^ ^^^i^'^ (pt) by fitting the Monte Carlo spectra to the data. 
The Monte Carlo spectra were referred to as total Monte Carlo spectra. The open 
charm contribution to a total Monte Carlo spectrum could be varied by using open 
charm differential cross sections parameterized by functions of the hadrons pr- The 
open charm cross sections were fit to the data using a least-squares minimization 
routine where each data point contributed an individual calculated as the square 
of the difference between the data and the total Monte Carlo spectrum for that data 
point, divided by the sum of the squares of the errors: 

,2 



X. - ^ , ^2 ' (5-1) 

where Nj is the number of muons from the data in a histogram bin j and W^'" is 
the number of muons determined from weighted spectra in bin j of a histogram of 
the total Monte Carlo spectrum corresponding to the data. The errors e and e are 
the errors for the data and composite Monte Carlo, respectively. 
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Use of single muons prevented the analysis from determining which species of 
hadron decayed to a given muon. In the case of open charm production, hadrons 
containing a charm quark and a light anti-quark (c, q) are designated D , and those 
having an anti-charm quark and a light quark (g, c) are designated by D. The inclusive 
D cross sections were extracted using /i+ spectra, and the inclusive D cross sections 
were extracted using /i~ spectra. 

Previous experiments have fitted various functions to their results based on 
theoretical predictions of the shape of the spectra as a function of px {h) } From the 
literature, the functions used to parameterize the open charm cross sections for this 
analysis were: 

Aiexp{-BpT) (5.2) 



(5.3) 



beam I 



jl-pT/Pbeamj ^^^^^ 

{p\ + a m^)" 

Another exponential form, exp {— B' p^) , used in |H] and P, was also at- 
tempted but resulted in very high total indicating it should not be used to describe 
the open charm cross sections for this data. Function 15.21 is referred to as the expo- 
nential form and the functions 15.31 and 15.41 are referred to as the 3 and 4 parameter 
forms, respectively. There were 4 fits performed using ^ ^ as the open charm 

cross sections, one where all parameters were free parameters, and one each where n 
was held fixed at the integer values 4, 5 and 6. These were chosen from theoretical 
predictions as well as fits used in previous experiments [5], EI], PHj - 

. Fits to Ao ^^V^J ^ ^^^'ITh were done for the same fixed values of n only. 



^References 0] - Jl] are a few examples. 
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All fits using the 3 and 4 parameter functions used = 1.5 GeV c ^. The 4 param- 
eter function includes an energy dependency term, -^^ere Pbeam is the center of 

Pbeam 

mass momentum of the proton beam. For this experiment, pbeam = 19.36 GeV/c. 

Each target provided four independent but related muon spectra; /i+ and 
fi~ from interactions in the target, and fi^ and /i~ from interactions in the dump. All 
four spectra were related by the number of protons incident on the target. Each fit 
required four total Monte Carlo spectra to be compared against the four data spectra 
simultaneously. 

A total Monte Carlo spectrum is the number of muons resulting from 2 x 10^ 
proton interactions in either copper or beryllium. More precisely, the total Monte 
Carlo copper target /i"*" spectrum is the number of fi~^ resulting from the interactions 
of 2 X 10^ protons in the copper target. The spectrum was determined from adding 
the number of fi^ from light hadron production, called the light contribution, plus 
the number of /i"*" resulting from production of open charm, called the open charm 
contribution. The light contribution is the sum of the muons from the production 
of 7r+, K+ and K^K Kl fi^ copper target Monte Carlo spectra. The light 
contributions were normalized since all the light hadrons were thrown using fully 
normalized differential cross sections. The open charm contribution was calculated 
using a set of values for the free parameters of the function being fitted. 

The four data spectra were related by the total number of incident protons, 
Np. The Monte Carlo generated all spectra using 2 x 10'' incident protons, so a 
scaling factor was introduced to scale the total MC spectra to the proper number 
of incident protons. The spectrometer lacked a beam intensity monitor which was 
sensitive to the requested number of protons per spill, so the parameter A^ (referred 
to as the scale factor A^) was introduced as an extra free parameter for all fits. 
The relationship between A^ and Ap is given in section 15.21 The requirement that 
the four data and four total MC spectra be used in calculating x^, and the use 
of normalized light contributions resulted in normalized open charm cross sections. 
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Figure 14.21 (page shows a representation of the process used to calculate the light 
MC contributions and the open charm transformation arrays used to extract the open 
charm cross sections. The open charm arrays and their use are described in Chapter 
15.31 For reference, Figure ITT] f Chapter IHl page 1 ^ shows the method used to analyze 
the single muon data and place the resulting muon spectra into arrays used in the 
least-squares minimization routines. 

Figure 15.11 shows the /i"*" spectra for the individual and the light Monte Carlo 
contribution used for all fits for the copper target (top left) and dump (lower left). 
The upper and lower right histograms show the light Monte Carlo contribution scaled 
to the data after the 3 parameter function with n a free parameter had been fitted. 
Figure EiH shows the /i"*" data (black circles), total Monte Carlo yU"*" spectra (blue open 
squares) and the open charm Monte Carlo contribution (red open circles) for the 
open charm cross section determined from the fit in Figure 15.11 Left is the spectra 
for the copper target and right is the dump. 



5.1 Minimization 

The three parameterizations of the open charm differential cross sections used in this 
analysis were: 

Ai exp {—Bpt) 

^ {I-PT / Pbeam)"' 

^ (Px + Oiml)^ 

The analysis histogrammed the data into 1-dimensional histograms with 40 
bins from to 10 GeV/ c in p^. The contents of the bins in the histograms were placed 
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Figure 5.1: Left Figures: light Monte Carlo fi^ contribution (black filled circles) 
from the copper target (top) and beam dump (bottom), composed of the sum of 
fi~^ spectra from vr"*" (red open circles), (blue open squares) and / (green 
open triangles). Right Figures: the same light Monte Carlo contribution after fitting 
the 3 variable function with N a free parameter (red open circles) and the /i"*" data 
(black closed circles). Errors are statistical only for the figures on the left. For 
figures on the right, the errors for the data are statistical only, the scaled light MC 
contributions have all systematic errors except normalization added in quadrature to 
the statistical errors. 
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Figure 5.2: Total open charm Monte Carlo /i"*" contribution (blue open squares), total 
light MC /i^ spectra from Figure 15.11 (red open circles) and the the data (black filled 
circles) from the copper target (left) and the dump (right). Data is the /i"*" spectra 
from the copper target, and the Monte Carlo spectra are from fitting the 3 variable 
function with n a free parameter. Data errors are statistical only. MC spectra have all 
systematic errors except normalization added in quadrature to the statistical error. 
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into arrays for use in minimization. The light Monte Carlo spectra were placed into 
similar arrays, and the contents of px if^) versus px (h) histograms were placed into 40 
X 40 2-dimensional arrays, referred to as hadron to muon pt transformation arrays. 

A change in any free parameter resulted in four total MC spectra to compare 
against the data. A ? called an individual ? x|) was calculated for all data 
points in the four data spectra using Equation 15. II fpage 17 ^) . The total ? xLt? 
the sum of the individual '■ 

k 

xlt = Y.^] (5.5) 

where k is the total number of elements in the four data arrays containing three or 
more muons.^ 

The spectrometer had the middle half of all x measuring hodoscope layers 
pulled, creating two inner acceptance edges (see Figure EiH page 1221) • These accep- 
tance edges combined with the single bend approximation used to correct for multiple 
scattering resulted in large uncertainties in muon px below 2.25 GeV/c. Attempts 
to correct this uncertainty failed because no decay vertex was available. As a conse- 
quence, the fits were limited to px {jj) > 2.25 GeV/c. 

The minimization routines used a grid search method to insure that local 
minima were avoided. The routines fitted each cross section individually in a ro- 
tating sequence. Normal initialization began with fitting the hadron cross section, 
E ^^^^ (Pt), while the anti-hadron cross section, E ^^^^ (Pt), remained constant 
using a set of initial parameter values. Once a set number of passes through the grid 
search had been performed on the hadron cross section, the routine re-initialized to 
begin fitting the anti-hadron cross section holding the hadron cross section constant 
at the values of the parameters that returned the smallest total (referred to as 
Xmin) during its fit. The scale factor N was allowed to float freely during fits to either 
cross section, with the requirement that all four total Monte Carlo spectra be scaled 

^The number of events in a bin needed to make a data point statistically significant varies between 
authors. This analysis chose the number of muons necessary to be 3. 
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by the same N simultaneously at each point in the routine where Xtot calculated. 

Because the target and dump were the same material in the copper target 
data, fits to the p-Cu data were performed first. Once the p-Cu data had been fit 
to the various functions, and the xLin values of the parameters of the p-Cu cross 
sections had been determined, fits of the p-Be cross sections were performed to the 
two p-Be target spectra. To accomplish this, the dump data were compared to scaled 
total Monte Carlo spectra calculated using a previous p-Cu result. The open charm 
contributions from the dump for fits to the p-Be data were calculated using the 4 
parameter function with n = 6 and the Xmin values of the parameters found from the 
fit to the p-Cu data. The scahng factor N remained a free parameter. 

5.2 Scaling Monte Carlo Spectra To The Target 
Material And Thickness 

The Monte Carlo spectra needed to be scaled by the integrated luminosity L, before 
they could be compared against the data, where 

^tgt = ^^^^ Xj [1 - exp {-I p I A/)] (cm"^) 

for the target spectra, and 

Ldmp = ^^'^^ \j exp {-Ip/ \i) (cm"2) 

for the dump spectra where Np is the total number of protons incident on the target, 
A^^ (mole~^) is the Avogadro constant. A/ (gm cm~^) is the nuclear interaction length 
for protons in the material, A (gm mole~^) is the atomic weight of the material, p (gm 
cm~^) is the density of the material, and / (cm) is the thickness of the target. The 
difference between the two is simply stating that all protons incident on the dump 
interacted in the dump, and the number of protons incident on the dump was the 
number incident on the target times the probability that the proton did not interact 
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while traversing the target material. The target scaling includes loss of protons due 
to interactions in the target material. 

It is important to note that Np is commmon to both. Since the experiment 
lacked a beam intensity monitor, the total number of protons incident on the target 
was unknown. To overcome this, the experiment split the total number of incident 
protons from the remaining terms of the two integrated luminosities such that L = 
Np L where L is 

Ltgt P I ^^y^ (^"^~^) 

for the target spectra, and 

^dmv ^ ^ e-xp {-I p / Xi) (cm"^) 

for the dump spectra. 

The minimization routine scaled the Monte Carlo spectra by the appropriate 
value of L to match the data being fitted at initialization. To estimate the number of 
incident protons Np, the total nuclear cross section weight thrown for p-Cu tt"*" (either 
the target or dump) is divided by the total p-Cu tt"*" cross section. The minimization 
routines input N as, & free parameter, and using the conversion from the thrown 
weight, Np ~ 21.42 x 10^1 N. The Xmin value for the parameter N thus determined 
the integrated beam flux corrected for dead time. 

5.3 Calculating The Open Charm Contributions 

The open charm contribution to a total MC spectrum was calculated using an array 
whose elements, W^~*'*, were the weight of the muons having transverse momentum 
between 0.250j and 0.250(j + 1) GeV/c, from open charm hadrons thrown with 

transverse momentum between 0.250i and 0.250(i + 1) GeV/c. The array used for 
calculating the /x"*" spectrum from and D*^ produced in the copper target was 
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composed of the elements 

Wg-^"(r, Cu) = Wf -^"(T, Cu) + Wf -^"(T, Cu) 

The Monte Carlo was used to develop a set of 6 arrays used to calculate the 
open charm contributions, one each fi~^ and fi~ from the copper target, the beryllium 
target and the copper beam dump. The dump required only one set of arrays because 
the difference between the arrays from the dump with either target presented to the 
beam would be the target thickness normalization presented in Chapter 15.21 The 
elements of the arrays were the values of the weight in the bins of 2-dimensional 
Pt (yu) versus Pt (h) histograms of the individual open charm hadrons thrown in 
the Monte Carlo. Figure 15.31 shows the 2-dimension p^ (yu) versus px {h) histogram 
resulting from 20 million and thrown from the copper target. Table I^TD gives 
the 12 individual arrays and the six transformation arrays used for calculating the 
open charm contributions. 



Z4/07/25 14.20 




Figure 5.3: The pr (/i) versus Pt{D) histogram of Monte Carlo D ^ used in the 
calculation of the muon contribution from D production from the copper target in 
the fitting routine. Weight does not include any cross section weight. 
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Table 5.1: The six transformation arrays used in the minimization routines to de- 
termine the open charm contributions from the copper target {Cu), beryllium target 
(Be) and dump (dmp). The contribution is determined from the transformation array 
which is the sum of the weights of the two individual arrays listed, after an arbitrary 
open charm cross section is applied. A fit required four contributions to be determined 
for each change in any of the free parameters. 



Spectrum 


Transformation 


Individual 


Desired 


Array 


Array 


fi+{Cu) 






H-{Cu) 


W'^^i'~{Cu) 


W^~^i'~{Cu) 


/x+(i?e) 






II- (Be) 






{dmp) 


W^^^'^{dmp) 




li~ {dmp) 


W^^^" {dmp) 


W^~-^^'~{dmp) 
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The value of an open charm (or anti-charm) cross section, W^^, was calculated 
at the center of each bin using p^^i = 0.250 (i — 1) + 0.125 (GeV/c) for i from 1 to 
40. The muon spectrum expected from throwing the open charm hadrons with the 
cross section being tested was then determined by 



40 



i=l 



Figure 15.41 shows the result of applying a cross section weight for to the 
histogram shown in Figure EISl 



Z4/07/25 '4.23 











■ c 





Figure 5.4: Same hadron to muon px distribution as in Figure loTal except trial input 
cross section weights W[^^^ have been calculated and applied to all columns i in 
Pt (h) during a minimization. The muon distribution Wj" expected from the input 
cross section is the projection onto the px (/i) axis. 



Figure 15.51 shows the projections of the histograms onto their respective axes 
to show the effect of applying a cross section (in this case an arbitrary D cross section 
found during a minimization). 
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Figure 5.5: The D spectrum of Figures 15.31 (top left) and 15.41 (top right) and the 
/i"*" spectra corresponding to the same Figures on bottom. All spectra are the projec- 
tion of the 2-dimensional pt (/^) versus pr [h) histograms onto their respective axes. 
Errors may be smaller than the symbols used, and include all systematic errors as- 
sociated with the decay and branching fraction parameterizations used to calculate 
those "weights as described in 15.41 
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5.4 Errors Used For Calculating ^ 



Only the statistical error was used for the data, 




Sj, where Sj is the total number 



of muons in element j. Errors for Monte Carlo spectra included the systematic and 
statistical errors of the parameterizations used to calculate the primary weights Wx, 
Wa, Wd, Wna and Wbr for light hadrons and Wd, Wna and Wbr for the open charm 
hadrons. See table 14. l| page ESI 

The total weight given an event in the Monte Carlo was defined as 



for open charm hadrons. Hadrons thrown from the dump had Wna = 1- The 
statistical error of s weighted events in a bin j is determined from |42j : 



All statistical or systematic errors of the parameterizations used in calculating 
the primary weights, ax, (Ja, ctd, ctna were also calculated. As an example, the 
systematic error, ax, of the differential cross section weight, Wx, of a light hadron 
was calculated as two parts, the errors of the parameters used to calculate the cross 
section, shown in Table fpage lU^. plus the normalization uncertainty, both of 
which were taken from the reference. The error due to the uncertainties in the values 
of the parameters was calculated by use of the error propagation equation, excluding 
covariant terms 



the light hadron differential cross section, Ai, B and m are three of those parameters 



Wk = Wx,kWAWDWNAWbr 



for light hadrons and 



Wk = Wd Wna W^r 



s 
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and is the error of the parameter Ai taken from the reference. The second error 
associated with the hght hadron differential cross section weight was the uncertainty 
in the absolute normalization, given by the reference as a relative error, R^. The 
error due to the normalization uncertainty was calculated by ax,n = WxRx- The 
error associated with the light hadron cross section weight, was then calculated 
as al = alp + a^ ,^. Additional errors for all other weights were calculated using the 
error propagation equation and the errors of variables from the references. The error 
on Wbr, (Tbr (strange and open charm hadrons only) was 0.2 percent (strange hadrons) 
and 11.5 percent (open charm hadrons) of the value of the branching fraction weight. 

Each muon from the Monte Carlo contributed the square of its total weight as 
a statistical error, plus the square of the error of the weights for use in calculating 



the e] 



where Wk is the total weight of muon k and a\ is the square of the errors of the 
weights 



a 



x,k 



+ 



a 



A,k 



a 



+ 



D,k _|_ 0" 



NA.k 

WWk 



+ 



^br,k 

Wlk 



The squared error, for s muons in a bin j is then 



fc=i 




(5.6) 



5.5 Parameter Errors 



Error routines were developed based on MINUIT, the method used for error analysis of 
function minimizations used by CERN . The errors of an n dimensional function, 
including all correlations, can be determined using contours of equal likelihood. By 
choosing the appropriate amount to add to the value of Xmin^ referred to as UP, 
the errors on the Xmin values of the parameters can be determined by minimization 
to n — 1 dimensions of the original function, where one parameter is held fixed. In 
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general, the method involves increasing or decreasing the value of the fixed parameter 
and then minimizing the remaining parameters until the smallest value of the total 
returned (referred to as Xtest here) remains smaller than or equal to xLax^ where 
Xmax = iXmin + UP). The maximum and minimum values the free parameters had 
to take in order that Xtest remain smaller than or equal to Xmax during the entire 
process were returned as their errors. 

The routines systematically increased/decreased the fixed parameters value by 
a pre-set amount until the returned value of Xtest increased to more than Xmax- Once 
this condition was met, the routine decreased the size of increment that the fixed 
parameter would increase/decrease and began minimization of the remaining free 
parameters where the value for the fixed parameter was the previous maximum or 
minimum value plus or minus the new increment size. This procedure was repeated 
until the increment size became smaller than 10"^. This experiment used values for 
UP for a 70% confidence level for both the nv = 3 and nv = 4 functions as given in 
Table O 
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Table 5.2: Values of UP for given confidence level and number of free parameters n. 
Taken from [IH] . 



Number of 
Parameters 


Confidence level (probability contents desired 
inside hypercontour of = xLin + ^P) 


50% 


70% 


90% 


95% 


99% 


1 


0.46 


1.07 


2.70 


3.84 


6.63 


2 


1.39 


2.41 


4.61 


5.99 


9.21 


3 


2.37 


3.67 


6.25 


7.82 


11.36 


4 


3.36 


4.88 


7.78 


9.49 


13.28 


5 


4.35 


6.06 


9.24 


11.07 


15.09 


6 


5.35 


7.23 


10.65 


12.59 


16.81 


7 


6.35 


8.38 


12.02 


14.07 


18.49 


8 


7.34 


9.52 


13.36 


15.51 


20.09 


9 


8.34 


10.66 


14.68 


16.92 


21.67 


10 


9.34 


11.78 


15.99 


18.31 


23.21 


11 


10.34 


12.88 


17.29 


19.68 


24.71 
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5.6 Other Errors 



This analysis did not require correcting the data for losses from acceptances, and 
the large error introduced by the addition of the two inside acceptance edges was 
minimized by cutting the first 5 data points from all spectra. Corrections for data lost 
due to data acquisition limitations, referred to as signal busy losses, are contained in 
the parameter N. Monte Carlo studies determined the uncertainty in the transverse 
momentum of the muons was less than 5 percent for the entire region above 2.0 
GeV/c. This uncertainty has the effect of shifting the pt spectrum to slightly higher 
values due to the fact that more events will 'feed down' than will 'feed up' because 
the spectrum is (approximately) exponentially decreasing as a function of pt- The 
same effect is present in all Monte Carlo spectra as well, so the relative effect between 
the two is, for the most part, cancelled. 

Loss of data due to inefficiencies in the hodoscopes was modeled in the Monte 
Carlo based on studies by the previous dimuon analysis. For that data all hodoscopes 
had efficiencies of over 96 percent, and for this analysis only the x measuring ho- 
doscopes were needed. Similarly, loss from bad wires and noise in the wire planes was 
modeled in the Monte Carlo spectra based on studies from the same analysis. As a 
check, the data were plotted using the variable TAN6'a;. This variable is very sensitive 
to the opening angle the muon trajectory had in the x — z plane, and any significant 
loss of events from a defective hodoscope or phototube should be clearly discernable. 
No such defects in the spectra were found. (See Figure ISHU page IHU) 

Light hadrons have a small uncertainty at moderate to high transverse momen- 
tum associated with the parameterization used. The fits determined by the BSC are 
valid up to approximately 6 GeV/cin The light hadron distributions above 6 

GeV/c have an additional uncertainty estimated to be less than 10 percent. This un- 
certainty can be neglected since the light hadron contributions are less than 1 percent 
beyond 6 GeV/c (see figure IHUZI) • 

Production of hadrons containing open bottom were ignored when calculating 
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the total Monte Carlo spectra. Semi-leptonic decays of open bottom hadrons intro- 
duces a small contamination since the decay may give rise to both a /i+ and a /i~ (the 
semi-leptonic decay ji^u^D^ where the then decays semi-leptonically to a 

/i~ is one example). The branching fractions of these modes are small, typically less 
than 3 percent, so contamination from production of open bottom is negligible. 

Table 15.31 lists the estmated uncertainty of the fits due to the analyzed muon 
transverse momentum, hodoscope inefficiencies, light hadron normalization, choice of 
the transverse momentum of the open charm hadrons at the mid-point value of the 
bin and contamination of the spectra through semi-leptonic decay of open bottom 
production. 

Table 5.3: Estimated errors on final cross sections from various sources not included 
in the errors of the parameters. The Muon and average hadron pt are pt dependent, 
and the error shown is the maximum for these errors at any px- 



Source 


Estimated 
Uncertainty (%) 


Muon Pt 


7 


Average Hadron pt 


5 


Hodoscopes 


3 


Light Hadron Normalization 


< 1 


Contamination From Bottom 


< 1 


Total 


9.2 
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Chapter 6 



Results 



This analysis used single muon data to extract the open charm differential cross 
sections as a function of the hadron transverse momentum. Eight variations of three 
functions were fit to the data; an exponential, a three parameter polynomial with n 
a free parameter, and both the three and four parameter polynomials with n fixed at 
the integer values 4, 5 and 6 (see Chapter page ITUj) . 

Table 16.11 gives the total number of muons in each data spectrum, the total 
number of muons in the data and the number of data points that each spectrum 
provided. Figure 16.11 shows the pt (a*) distributions from the copper target (top) 
and beryllium target (bottom). The vertical line represents the minimum transverse 
momentum used (2.25 GeV/c) during fitting. 



6.1 Fits To The Data 

The values for the parameters for the best fits to the data using the exponential func- 
tion are presented in Table WJA for the 3 parameter polynomial function in Table 
and the four parameter polynomial in Tables IHTil (copper) and 16.51 (beryllium). The 
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Table 6.1: The number of muons, N^, and the number of data points, Ndp, after 
final cuts were applied for all spectra used in this analysis. T indicates events from 
the targets, and D indicates events from the dump. 



S]:)ectrum 






T+ 


16657 


14 


T- 


11954 


14 


Copper -D+ 


59576 


20 


D- 


53322 


18 


Total 


141689 


66 


r+ 


13368 


16 


T- 


9482 


14 


Beryllium D'^ 


120262 


19 


D- 


109017 


19 


Total 


252129 


68 
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Figure 6.1: The single muon E866 data spectra. Top two figures are the /i+ (red 
triangles) and fi~ (blue upside down triangles) for the copper target data set from 
the target (left) and dump (right). Bottom row is the spectra from the beryllium 
target data set. Vertical lines show lower limit (pt (a*) > 2.25GeV/c) of the muon 
transverse-momentum used during fitting. Both dump figures are the events from 
the copper dump when the designated target (Cu or Be) was presented to the beam. 
Errors are statistical only. 
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total x^, Xmm' ^^"^ reduced x^, ^prf/' shown for the fit for each parameteriza- 

tion. From the tables, the author concludes: 

• The 3 parameter function with either n a free parameter or n fixed at n = 6 as 
well as the 4 parameter function with n fixed at n = 6 provided the best fits to 
the data. 

• Though no errors were calculated for floating values of n, the weighted average 
value of n for all fits with n a free parameter was 6.18 ± 0.19. 

• The exponential function provided a remarkably good fit for a large range in 
PT{h),2.Q<pT <7.0GeV/c. 

• The 4 parameter function provides no additional information beyond that de- 
termined from fitting the 3 parameter function. 

All fits of functions having 4 free parameters, the 3 parameter floating n and all 
4 parameter functions, had large correlations between the parameters. No errors were 
calculated for fits using the 3 parameter function with n a free parameter. Because 
of the large correlations, two sets of errors were presented for fits using the fixed n 4 
parameter function; one where the parameter a was held constant at its Xmm value 
(referred to as the fixed a errors), and one set where the parameter m was held fixed 
at its xlain value (fixed m errors). 

Figures 16.21 through 16.91 show the total Monte Carlo spectra plotted against 
the data for four selected fits: the exponential function, the 3 parameter polynomial 
function with n free and n fixed to = 6, and the 4 parameter polynomial function 
with n fixed ton = Q. All Monte Carlo spectra are those at Xmm function being 

fitted. Errors for the Monte Carlo spectra include all statistical plus all systematical 
errors as explained previously. The Monte Carlo errors do not include the errors 
for the parameters of the input open charm corss sections. All data spectra have 
statistical errors only. Figure 16.101 is the same as figure 16.61 except the open charm 
and light contributions are included. 
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Table 6.2: Parameter values from fits to the data using Ai exp {—Bpt) as the input 
shape of the open charm differential cross section. Values for D are top line and D 
are the bottom line. Errors are the 70 percent confidence level errors as described in 
Chapter 15.51 



Copper Target Data 


Ai B N 


Arrwn 


Xpdf 


1 qi +0.32 . ni +0.04 

L.JL _Q27 ^-^^ -0.04 0.23 

1 SQ +0.34 . 04 +0.04 • 0.47 
x.oo _Q28 -L--''* -0.04 


81.4 


1.3 




Beryllium Target Data 


Ai B N 


A.min 


Xpdf 


0.60 +'f, 2.11 

-0.23 -0.12 0.10 
n +0.16 9 1 +0.07 ■ 0.66 
U.O^ _p21 ^-J-O -0.12 


62.9 


1.0 
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Table 6.3: Parameter values for D (top line) and D (second line) from fits to the data 
where the differential cross section was input as / „ ^"^ ^.n . Top two lines of both 

target data sets are the full float minimization. Errors are calculated for 70 percent 
confidence level as explained in Chapter 15.51 



Copper Target Data 

M a n N xldf 


649710 4.14 6.32 

21.84 

286500 3.77 6.19 


54.09 


0.92 


37.52 0.07 26.03 

— 3.45 —0.07 ^ Q —0.68 

31.91 !Jf 0.06 ■ 26.03 tjg 


305.41 


5.01 


2334 1.71 22.39 tj;? 

1931 1.66 tlf, ■ 22.39 


99.58 


1.63 


161765 ;S 3-53 tf^ 21.70 
130126 3.44 tin ' 21.70 


56.08 


0.92 




Beryllium Target Data 

^2 a n N xLn xldf 


19238 3.03 6.06 

42.19 

17013 2.92 6.14 


60.14 


0.99 


6.50 0.08 42.39 

4.88 tiil 0.08 tZ ■ 42.39 tlf. 


167.32 


2.66 


289 +S 1.46 41.97 +»2 

— 75 —0.39 _ „ —0.46 

99f- +13 . p-r, +0.08 ■ Q7 +0.24 
ZZO i.OU _Q.^3 ^i.y^ _Qg2 


84.44 


1.34 


15280 2.94 +2:^^ 42.13 tlf, 
11280 2.91 • 42.13 


62.19 


0.99 
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Table 6.4: The Xmin parameter values when the open charm differential cross section 
was input in the form and fitted to the copper target data. Top line is 

for D and second line is for D. All fits had n fixed to the values shown, pbeam = 19.38 
GeV/c for this analysis. Errors were calculated holding a fixed at its Xmin value (top 
pair of lines for each n) and holding m fixed at its Xmin value (bottom pair of lines for 
each n) because of strong correlations. Errors are calculated for 70 percent confidence 
level as explained in Chapter 15.51 



Copper Target Data 

^2 a n m N xLn xldf 


1412 +1^° 2 75 8 00 +°-^^ +°-^^ 

1616 ^lll 2.99 8.80 +Sg ^If^ 

4.0 ' 21.02 ~ 
1412 +158 2 75 8 00 

-1^-^^ -156 -^-'"^ -0.24 -0.23 

1616 2 99 +°-^^ 8 80 +°-^° 


58.48 


0.99 


1^9^9 +^i^° 3 30 4 17 

insn/i +1520 Q no q si +0-5i +0.10 

10804 _^283 '^•"^ "^-^l -0.51 -0.24 

ic;qc;q +1940 on +0.20 ..7 ^1.95 ^^^^ 

iOJOJ -^.^^Q 21 ^--L' -0.35 

10804 3.02 3.8I 


54.61 


0.93 


300220 4.07 1.29 

137891 -™ 3.48 0.13 

300220 4.07 " 1.29 ' 

137891 3.48 ^J;- 0.13 


54.23 


0.92 
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Table 6.5: The xLin parameter values when the open charm differential cross section 
was input in the form Ao and fitted to the beryllium target data. Top 

line is for D and second line is for D. All fits had n fixed to the values shown. 
Pbeam = 19.38 GcV/c for this analysis. Errors were calculated holding a fixed at its 
Xmin value (top pair of lines for each n) and holding m fixed at its Xmin value (bottom 
pair of lines for each n) because of strong correlations. Errors are calculated for 70 
percent confidence level as explained in Chapter 15.51 



Beryllium Target Data 

^2 a n m N xLn xldf 


508 +;g 3.34 11.23 

402 3 1^ 1187 
402 3.15 11.87 

508 +^^^ 3 34 +°-^^ 1123 +°-°^ 

/,p,r, +126 o 1 cr +0.60 07 +0.04 
"iUZ J.iO _Q5g il.O^ 


59.30 


0.97 


4151 +1228 04^5 f5 CTQ +1.15 +0.10 
^^^^ -903 "-""^ -1.06 -0.45 

q/|oq +1157 007 79c: +1.30 +0.10 
<3^^'3 _g^Q 0.0^ ^.ZO -^-^g 

4ic:i +1167 nAa +0.46 ^ c-q ^^.18 ^^^^ 
^^^^ -946 -0.45 "-""^ -0.45 
'^/IS'^ +1105 007 +0.50 79c- +0.05 


58.44 


0.96 


R7c;Q7 +19855 001 +1.13 +0.19 
-14577 ^--^-^ "^"^-^ -1.05 -0.45 

68581 +™ 4 19 4 49 +°-'' 

uouo± -16101 ^--^^ ^-^^ -1.20 .,,.0 
(^7f-q7 +20187 A .1 +0.43 001 42. iS +0.10 
KiiO^l _i6248 -0.43 '^''^^ -0.67 

68581 +2^^^^ 4 19 4 49 +°-°^ 

UOJOI _i8335 _o.48 ^-^^ -0.44 


58.25 


0.95 
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Figure 6.2: The copper target data (black closed circles) and total Monte Carlo 
spectra at Xmin (^^^ open circles) from the fit where the D/D cross sections were 
of the form Aiexp{—Bpx)- Histograms are /x"*" (left) and n~ (right) for the copper 
target (top) and copper dump (bottom). Errors for the data are statistical only. 
Errors for the Monte Carlo include all errors used for calculating x^(see text). Error 
bars may be smaller than symbols used. 
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Figure 6.3: The beryllium target data (black closed circles) and total Monte Carlo 
spectra at Xmin i^^^ open circles) from the fit where the D/D cross sections were of 
the form Aiexp{—B pj). Histograms are /i+ (left) and ^~ (right) for the beryllium 
target (top) and copper dump (bottom). Errors for the data are statistical only. 
Errors for the Monte Carlo include all errors used for calculating x^(see text). Error 
bars may be smaller than symbols used. 
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Figure 6.4: The copper target data (black closed circles) and total Monte Carlo 
spectra at Xmm (^^d open circles) from the fit where the D/ D cross sections were of 
the form / „ ^xn . The exponent n was a free parameter for this fit. Histograms 

are /i+ (left) and fi~ (right) for the copper target (top) and copper dump (bottom). 
Errors for the data are statistical only. Errors for the Monte Carlo include all errors 
used for calculating x^(see text). Error bars may be smaller than symbols used. 
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Figure 6.5: The beryllium target data (black closed circles) and total Monte Carlo 
spectra at xLin (^^d open circles) from the fit where the D/D cross sections were of 
the form / „ ^xn . The exponent n was a free parameter for this fit. Histograms are 

fj,~^ (left) and fj,~ (right) for the beryllium target (top) and copper dump (bottom). 
Errors for the data are statistical only. Errors for the Monte Carlo include all errors 
used for calculating x^(see text). Error bars may be smaller than symbols used. 
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Figure 6.6: The copper target data (black closed circles) and total Monte Carlo 
spectra at Xmm (^^d open circles) from the fit where the D/D cross sections were 
of the form / „ ^xn . The exponent n was fixed at n = 6 for this fit. Histograms 

are /i"*" (left) and /i~ (right) for the copper target (top) and copper dump (bottom). 
Errors for the data are statistical only. Errors for the Monte Carlo include all errors 
used for calculating x^(see text). Error bars may be smaller than symbols used. 
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Figure 6.7: The beryllium target data (black closed circles) and total Monte Carlo 
spectra at xLin (^^d open circles) from the fit where the D/D cross sections were of 
the form / „ ^-.n . The exponent n was fixed at n = 6 for this fit. Histograms are 

/i"*" (left) and fj,~ (right) for the beryllium target (top) and copper dump (bottom). 
Errors for the data are statistical only. Errors for the Monte Carlo include all errors 
used for calculating x^(see text). Error bars may be smaller than symbols used. 
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Figure 6.8: The copper target data (black closed circles) and total Monte Carlo 
spectra at Xmm (^^d open circles) from the fit where the D/ D cross sections were of 
the form A2 (^~Pjr/pfteam) ^ exponent n was fixed at n = 6 for this fit. Histograms 

are /i+ (left) and /i~ (right) for the copper target (top) and copper dump (bottom). 
Errors for the data are statistical only. Errors for the Monte Carlo include all errors 
used for calculating x^(see text). Error bars may be smaller than symbols used. 
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Figure 6.9: The beryllium target data (black closed circles) and total Monte Carlo 
spectra at xLin (^^d open circles) from the fit where the D/D cross sections were of 
the form A2 ^^T'^T' ^ ^'""tT^ ■ The exponent n was fixed at n = 6 for this fit. Histograms 

are /i+ (left) and /i~ (right) for the beryllium target (top) and copper dump (bottom). 
Errors for the data are statistical only. Errors for the Monte Carlo include all errors 
used for calculating x^(see text). Error bars may be smaller than symbols used. 
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Figure 6.10: The copper target data (black closed circles) and total Monte Carlo 
spectra at Xmin (^^^ open circles) from the fit where the D/D cross sections were of 
the form . „ ^.n . Blue triangles are the open charm contribution and the green 

(p2+am2) to 1- to 

upside-down triagles are the contribution from light hadrons. The exponent n was 
fixed at n = 6 for this fit. Histograms are /i"*" (left) and fi~ (right) for the copper target 
(top) and copper dump (bottom). Errors for the data are statistical only. Errors for 
the Monte Carlo include all errors used for calculating x^(see text). Contributions 
contain all systematic errors used in calculation of the weights added in quadrature 
to the statistical error, but do not include errors of the fit. Error bars may be smaller 
than symbols used. 
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6.2 The Open Charm And Anti-charm Cross Sec- 
tions 



The differential p-Cu and p-Be open charm cross sections, determined from fits to the 
data shown in figures through 16 .91 are shown in figures lU. Ill and f6. 121 Errors were 
calculated using the error propagation equation as described in 15 .41 for the parameters 
in the function used to describe the cross section added in quadrature to the error 
of the scaling variable N . The error shown represents a 70 percent confidence level 
that if the value of the first parameter was held within the stated range in values, 
the cross section would lie between the upper and lower error bars. Errors for the 
4 parameter function were performed twice because of large correlations; one set of 
errors were calculated holding the value of a fixed (solid lines), and one set with m 
fixed (dashed lines). Fits to the 3 parameter function with n a free parameter had 
large correlations as well, and no errors were calculated for those fits. For clarity, the 
same cross sections are plotted without errors in figure The p-Cu and p-Be D 
and D differential cross sections are plotted together for the four selected functions 
in figures I(i.l4| and 16.151 The four fits selected appear to be well constrained, even 
when projected to = GeV/c, as shown in figure IHTTHl Figure IU'.17I shows the 
D differential cross section scaled hy A = 63.546 plotted with the vr"*" and cross 
sections given in j22]- 
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Figure 6.11: The differential cross sections for D (left) and D (right) determined from 
fits to the copper data for the various functions: exponential (black), 3 parameter 
with n a free parameter (red), the 3 parameter with n fixed to n = 6 (blue) and 
the 4 parameter with n = 6 (green). Errors shown include the error on added in 
quadrature. 
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Figure 6.12: The differential cross sections for D (left) and D (right) determined from 
fits to the beryllium data for the various function: exponential (black), 3 parameter 
with n a free parameter (red), the 3 parameter with n fixed to n = 6 (blue) and 
the 4 parameter with n = 6 (green). Errors shown include the error on N added in 
quadrature. 
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Figure 6.13: The p-Cu (top) and p-Be (bottom) D (left) and D (right) differential 
cross sections from fits to the data for the various functions: exponential (black), 3 
parameter with n a free parameter (red), 3 parameter with n = Q (blue) and the 4 
parameter with n = Q (green). No errors are shown for clarity. 
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Figure 6.14: The p-Cu D (black) and D (red) differential cross sections for the various 
functions: exponential (top left), 3 parameter with n a free parameter (top right), 3 
parameter with n fixed at n = 6 (bottom left) and the 4 parameter function with n 
fixed at n = 6. Errors include the error on the parameter added in quadrature. 
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Figure 6.15: The p-Be D (black) and D (red) differential cross sections for the various 
functions: exponential (top left), 3 parameter with n a free parameter (top right), 3 
parameter with n fixed at n = 6 (bottom left) and the 4 parameter function with n 
fixed at n = 6. Errors include the error on the parameter N added in quadrature. 
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Figure 6.16: The p-Cu D differential cross section projected to low pt for the various 
functions: exponential (black), 3 parameter with n a free parameter (red), 3 parameter 
with n fixed at n = 6 (blue) and the 4 parameter function with n fixed at n = 6 
(green). No errors are shown. 
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Figure 6.17: The p-Cu D differential cross section projected to low pt for the 3 
parameter function with n = 6, and scaled down hj A = 63.546 (blue). Black is the 
TT"*" and red is the differential p-p cross sections as given in |32]. 
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6.2.1 Comparison To Other Experiments 

The open charm cross sections from this analysis were extracted from single muon 
events from 2.25 to ~7 GeV/c. Extensive Monte Carlo studies have determined that 
contributions from open charm hadrons to the single muon spectra had 0.2 ^ Xp ^ 0.8 
and 2.25 ^ Pt ^ 7.0 (GeV/c). The open charm cross sections were extracted using 
3 functions to describe the shape of the cross section: Ai exp {—Bpt), / ^ ^"^ ^yu and 
A2 , Comparison to other p-p or p-A results at similar energies is prob- 

lematical since all fit their data to functions dissimilar to those chosen here. The 
increased range in transverse-momentum from this analysis is far beyond other p-A 
experiments at this energy, covering a much larger range in as well. 

The LEBC-MPS Collaboration (E743) [6 fitted their 800 GeV/c p-p data 
with the function 

and E789 |^ fitted their results from 800 GeV/c p-A data with the function 

Pt c'^'Pt 

The exponential using the square of the transverse-momentum failed to adequately 
describe the data from this analysis because it falls too steeply over the range in 
transverse- momentum. 

Figure 11.51 shows the per nucleon total inclusive neutral open charm cross 
section from E789, and Figure 11.41 shows the p-p total inclusive open charm cross 
section from E743. For comparison to the results from this analysis, the E743 data 
was divided by 2 since their results were the total open charm cross section which 
is assumed to be twice the neutral cross sections determined by E789 and twice the 
sum of the D and D cross sections found by this analysis, and both the E743 and 
E789 results were scaled by Aqu = 63.546 and Ase = 9.012182. The sum of the D 
and D cross sections determined from fitting the 3 parameter function with n = 6 to 
the p-Cu data is plotted with the scaled results from E743 (left) and E789 (right) 
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in Figure 16.181 and the result from fitting the same function to the p-Be data is 
shown in Figure 16.191 The reader is cautioned that the cross sections shown from this 
analysis below 2.25 GeV/c are the projected cross sections, since the analysis used no 
data below that momentum. The reader is also cautioned that both the E743 and 
E789 data are scaled to A assuming a^px) = 1.0 (this is a result presented in E789). 
These comparisons have also ignored the ratio of charged to neutral production. These 
figures show the results from this analysis lie between the two previous measurements. 
This confirms that the use of single muon spectra and the methods adopted by this 
analysis to extract the differential cross sections provided reasonable results. 

Table inini provides results from selected experiments. It lists the type of data, 
Pt and xp ranges covered, the function(s) and values of the parameters from fits to 
their data and results from this analysis. From figures 16.181 16.191 and the table the 
author concludes: 

• The open charm cross sections from this analysis reproduce the open charm 
Pt distributions from both previous 800 GeV/c p-p and p-A experiments, 
when projected to low hadron p^- 

• Data from 250 GeV/c p-A interactions (E769 [9 ) was fit to the exponential 
function. The value of B reported was 3.0 ± 0.3, while the value of B found 
by this analysis was 1.91 ± 0.04 (copper) and 2.11 ± 0.12 (beryllium). This is 
reasonable since the E769 data was from 0.0 < pr ^ V^O GeV/c, while the 
data from this analysis is at higher p^. The slope of the cross section becomes 
smaller with increasing pt , as evidenced by the results from this experiment. 

• 250 GeV/c n-A data from E769 was fitted with the 3 parameter function, where 
a = 1.4±0.3 and n = 5.0±0.6. This analysis determined that n = 6 and a ~ 3.5 
(copper) and a ~ 2.9 (beryllium), indicating that open charm production, as 
a function of Pt, is softer for proton interactions than meson interactions (see 
figure II. 6|) . 
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Table 6.6: Top is the experiment, interaction, beam momentum and range in pt and 
xp covered. Bottom table gives the parameterizations used to fit the data by the 
experiment and the values of the parameters reported. Note that the value of B' from 
E789 does not include all pt dependencies of the function used by that experiment 
(see text). Values given for this analysis (E866) are for fits to determine the D cross 
section only. N/A indicates that the parameterization was attempted but resulted in 
very large Xmin values of the parameters are not shown. Comparison of the 

cross sections from this analysis and the E743 and E789 experiments are shown in 
figures 16.181 and 16.191 
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Figure 6.18: Comparison of the results from this analysis from fitting the p-Cu data 
using the 3 parameter function with n = 6 to the E743[n] (left) and E789[S] (right) 
data. The cross section from this analysis is the sum of the D and D cross sections. 
The invariant cross sections for E743 and E789 were determined by dividing their 
published results by (1 — |xi?|)^'^ x 2 x tt. The data were then scaled by Acu = 
63.546. Black vertical line indicates the minimium pt (/i) data used by this analysis. 
The reader is cautioned that the data from E743 and E789 are evaluated at xp = y = 
and compared to the cross section from this analysis projected to lower momentum. 
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6.2.2 The Ratio Of Charged To Neutral Production 

Theory suggests that 

where Rcn is the ratio of charged to neutral open charm production j3] ^2]- Prior to 
2000, experimental results from proton induced charm production showed the ratio 
to be consistent with unity [H] (Zj. Results from the production of open charm via 
pion-nucleon interactions show the ratio from 0.27 to 0.50.^ The discrepancy between 
the ratios from meson interactions versus proton interactions, which should be the 
same, has generated some interest in the community A more recent study by the 
HERA-B collaboration [13] found the ratio to be 0.54 ± 0.11 ± 0.14 in 920 GeV/c 
p-A collisions. Table 16.71 shows four experimentally determined charged to neutral 
production ratios for open charm from p-p and p-A interactions. 



Table 6.7: The charged to neutral production ratio, Rcn, from other p-p or p-A 
experiments. 



Experiment 




Momentum 
GeV/c 


Ratio 


E783 ^ 


p-p 


800 


1.2 ± 0.6 


E653 


p-A 


800 


1.0 ± 0.6 


LEBC-EHS HE 


p-p 


400 


0.7 ± 0.1 


HERA-B [in! 


p-A 


920 


0.54 ± 0.18 


Weighted Avera 


ge 


0.7 ± 0.1 



^ These results are presented in tabular form in [S] 
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Beryllium Target 

E743 E789 




2468 10 02468 10 

Pt (GeV/c) 

Figure 6.19: Comparison of the results from this analysis from fitting the p-Be data 
using the 3 parameter function with n = 6 to the E743[n] (left) and E789[S] (right) 
data. The cross section from this analysis is the sum of the D and D cross sections. 
The invariant cross sections for E743 and E789 were determined by dividing their 
published results by (1 — I^f])^'^ x 2 x tt. The data were then scaled by Ase = 
9.012182. Black vertical line indicates the minimium pt (h) data used by this analysis. 
The reader is cautioned that the data from E743 and E789 are evaluated at = ?/ = 
and compared to the cross section from this analysis projected to lower momentum. 
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This analysis lacks the sensitivity required to measure the charged to neutral 
production ratio, Rcn- This is due to the large transverse momentum shift from 
the parent hadron to the decay muon. This shift resulted in the shapes of the muon 
spectra being very similar for all open charm hadrons. During fitting the contributions 
from both the charged and neutral open charm mesons were added together, to give 
the total open charm Monte Carlo contribution. The analysis used Rcn — 1 for all 



E ^ (pt) and E '^^^^ (Pt)- The D cross section was then the cross section for 



fits, so the two open charm cross sections determined from each fit were designated as 

either or D^, and the D cross section was the cross section for either the D~ or 

By using a charged to neutral production ratio other than 1, fitting would then 
be to either the charged cross section, or the neutral cross section, depending on how 
the ratio was introduced into calculating the open charm and open anti-charm Monte 
Carlo spectra. An open charm Monte Carlo spectrum is used as an example. The 
open charm contribution was determined from 

40 
i=l 

Choosing aD+ — Rctdo the equation above may be re- written as 

40 

i=l 

Further, assuming the two contributions are related by (for this anlysis, the open 
charm contributions had virtually identical shapes, though different magnitudes) 

where s is a constant, the open charm contribution to the total Monte Carlo spectra 
can be expressed as 

40 



W^^^{sR+l)J2 [Kf "^fj^'"] (6-1) 

i=l 
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While it is possible to include the ratio as a free parameter in the fitting routines 
in much the same manner as the free parameter N, the similarity in the shapes of 
the muon spectra from D'^ and production make the sensitivity to Rcn very low. 
The free parameter N is constrained by the light contribution which is different in 
shape to the open charm contribution. 

Effects of Rcn 7^ 1 can be examined by using Equation 16.11 above. Values 
of -Rc77vless than one result in a smaller open charm contribution, which results in 
a cross section larger than the one determined with Rcn = 1, and values of 
i^c* AT greater than one result in cross sections that are smaller. The amount they 
would differ is given by 

Figure lU. 201 shows the D cross section from the fit to the p-Cu data where the 
shape of the cross section was assumed to be given by the 4 parameter function with 
n fixed at = 6 (red). The errors shown are those calculated holding a fixed at its 
Xmin value. The blue line is the cross section, and the green line is the cross 
section that would be found for Rcn = 0.7. 
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Difference Between R^^ = 1 And Rq^ = 0-7 



.16-1 



.16-2 



> -16-3 
o 

n 
U 

El 6-05 



CO 

Q. 

b 1 6-06 

CO 

■D 
lU 

16-07 



16-08 



16-09 




5 6 
Pt (GeV/c) 



10 



Figure 6.20: The D cross section determined by fitting the 4 variable function with 
n fixed at n = 6 and Rcn = 1 to the p-Cu data (red), and the D° (blue) and 
D+ (green) cross sections that would result for Rcn = 0.7. All errors are calculated 
with the value of a held fixed to its Xmin value. 
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6.2.3 Nuclear Dependency 

Nuclear dependency is defined by 



aA = A"' aN 

where a is a function of either xp, Pt ot both, A is the atomic weight of the material 
and N is nucleon. E789 studied the production of neutral open charm produced in 
800 GeV/c p-A collisions near xp = 0. The nuclear dependency reported by this 
study was a = 1.02 ± 0.03 ± 0.02, implying that the cross sections scaled as the 
number of nucleons In contrast, E866 |23| measured the power a for hidden 
charm for three regions in xp, SXF (—0.1 < xp < 0.3), IXF(0.2 < xpO.6) and LXF 
(0.3 < xp < 0.93). For this analysis, 0.2 < a;^ < 0.8, which lies between the IXF and 
LXF data from E866. Their results for the power a{pT) are shown in Figure lU. 211 
The ratio _ 

_ ABe E J ' (pt) 

AcuE ^) ' (pt) 

for the 3 parameter n = Q results from this analysis are shown in Figure 16.221 The 
power a can be determined from the ratio of the cross sections, r = (ycu/<^Bei and the 
relations acu = and ase = "^-^v- 

This analysis has extracted a{pT) for both ratios {r{D) and r{D)), and the 
resulting a{pT) are shown with the IXF and LXF J/\I' data from [23J in Figure lU. 231 
Both ratios used the fits of the 3 parameter n = 6 to the copper and beryllium data. 
Based on these results, this analysis concludes that the cross sections, as a function 
of Pt, do not simply scale by the atomic weight A (admittedly within large errors, 
though). 
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0.7 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 
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Figure 6.21: a versus for J/^ (solid circles) and \E'' (open boxes) production by 
800 GeV/c protons. Results are shown for the three data sets - SXF, IXF, and LXF 
(see text) - which have {xp) = 0.055, 0.308 and 0.480 respectively. Only statistical 
uncertainties are shown. An additional systematic uncertainty of 0.5% is not included. 
Also shown are the NA3 results at 200 GeV/c [51J. The solid curves represent the 
parameterization a{pT) = Ai{l + 0.0604pr + 0.0107p|,), where Ai = 0.870, 0.840, 
0.782, and 0.881 for the SXF, IXF, LXF data sets, and the NA3 data, respectively. 
Taken from [SB] 
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Figure 6.22: The ratio of the copper and berylhum D cross sections (black), and the 
copper and berylhum D cross sections (red). Both ratios are normalized to Ase/^cu 
and use the results of the fits to the data of the 3 parameter polynomial function 
with n = 6. Thin black (red) lines are the scale of the errors of the ratios. Errors 
include the error of N added in quadrature to the statistical and systematic errors of 
the other parameters. 
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2^ 4 6 8 ^0 



Pt (GeV/c) 

Figure 6.23: The power a{pT) determined by this analysis for D + D, where the 
results from fits using the 3 parameter polynomial with n = 6 were used as the open 
charm/anti-charm copper and beryllium cross sections. Red open circles are the IXF 
data, and the blue open circles are the LXF data taken from [23]. Thin black lines 
are the errors for a from this analysis. Errors include the error for the parameter N 
added in quadrature to the error of the remaining parameters. 
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6.3 Summary 



The author takes the position that determination of the 'shape' of the open charm 
cross sections is more important than determination of the absolute magnitude. 
Clearly both would be optimal, but the use of secondary spectra required making 
several assumptions such as the ratio of charged to neutral open charm production. 
These assumptions provided uncertainties in the absolute magnitudes of the open 
charm and open anti-charm cross sections that could not be removed. 

Single muon data from 800 GeV/c p-Cu and p-Be interactions were used 
to extract the differential open charm/ ant i- charm cross sections as a function of px- 
Several functions describing theoretical predictions of the shapes of the cross sections 
were fit to the data. Monte Carlo studies indicated that the open charm contributions 
were from open charm hadrons produced with 2.25 < pr ^ 7.0 (GeV/c), < z/cm. ^ 
2.0 and 0.2 ^ xp ^ 0.8. Tables 16 . 31 and 16 . 51 show that the 3 parameter (with n 
a free parameter or fixed at the integer value 6) and 4 parameter (with the parameter 
n fixed at the integer value 6) polynomial functions resulted in the smallest Xpdf ^"^^ 
both the p-Cu and p-Be data. Extrapolation of the results from the 3 parameter 
function with n = 6 show good agreement with previous experiments. Fits of the 
three functions to both the p-Cu and p-Be data indicate that the parameter n is 
very close to 6.0. Due to large correlations introduced by the scaling factor and 
loss of data below 2.25 GeV/c, errors for the parameters from fits using the 3 variable 
function with n a free parameter cannot be determined accurately. In summary: 

• The production of open charm hadrons, as a function of the hadron p^, is well 
represented by the function: 

(py + ami)"' 

with n = 6. Previous fits using this function have only been reported for tt-A 
data, where it was found that the same function fit the data well, and n = 5.0. 
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• The cross sections can also be well described using the simple exponential 

Ai exp{-BpT) 

for 2.25 < pr < 8.0 GeV/c. 

• The weighted average value (over D and D) of a for the 3 parameter function 
and n = 6 for p-Cu production is 3.49 ± 0.01, while the weighted average for 
p-Be production is 2.92 ± 0.02. The difference indicates that the production 
via p-Cu has a steeper slope than p-Be production for the region in pt covered 
by this analysis. 

• Extraction of the power a used to relate proton-nucleon cross sections to proton- 
nucleus cross sections, often referred to as the nuclear dependency, indicates a 
structure similar to that found for the production of J/\E', though large uncer- 
tainties at low Pt make it difficult to claim that such a structure was found 
with certainty. Regardless, the author claims that, for the region in pt for this 
analysis, the power a continuously rises even if the large errors are taken into 
account at pt ~ 7.0 GeV/c. The conclusion reached for this analysis is that 
the cross sections do not simply scale by the nuclear weight A. 
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Appendix A 

Miscellaneous Apparatus 
Information 



Table A.l: Drift chamber information, z is given for the average of the two planes 
for each pair. Chambers are listed in ascending distance along the z axis. All lengths 
are in inches. 



Detector 


Angle 


Number 


Wire 


z 


Aperture 


Operating 


Pair 


(deg) 


of Wires 


Spacing 


{x X y) 


Voltage 


Vl-Vl' 


-14 


200 


0.250 


724.69 


48.0 X 50.0 


+1700 


Yl-Yl' 





160 


0.250 


748.81 


48.0 X 40.0 


+1700 


Ul-Ul' 


+14 


200 


0.250 


755.48 


48.0 X 50.0 


+1700 


U2-U2' 


+14 


160 


0.388 


1083.40 


66.0 X 62.1 


-1950 


Y2-Y2' 





128 


0.400 


1093.21 


66.0 X 51.2 


-2000 


V2-V2' 


-14 


160 


0.388 


1103.25 


66.0 X 62.1 


-2000 


U3-U3' 


+14 


144 


0.796 


1790.09 


106.0 X 114.7 


-2200 


Y3-Y3' 





112 


0.820 


1800.20 


106.0 X 91.8 


-2200 


V3-V3' 


-14 


144 


0.796 


1810.24 


106.0 X 114.7 


-2200 
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Table A. 2: Information on the banks of hodoscopes. Planes are presented in the order 
in which a particle would traverse them. Minimum and maximum values oi\P ^1 "P z\ 
for each x plane are presented for reference, and are based on disconnecting the high 
voltage supplies to the middle half of each bank as shown in figures EISl and All 
lengths are in inches. 







Number 


Scint. 


Aperture 




TAN^^. 




z 














Scint. 


Width 


X X y 


Min 


Max 


Yl 


769.78 


16 


2.5 


47.50 X 40.75 






XI 


770.72 


12 


4.0 


47.53 X 40.78 


0.015 


0.030 


Y2 


1114.94 


16 


3.0 


64.625 X 48.625 






X3 


1822.00 


12 


8.68 


105.18 X 92.0 


0.014 


0.030 


Y3 


1832.00 


13 


7.5 


104.00 X 92.0 






Y4 


2035.50 


14 


8.0 


116.00 X 100.00 






X4 


2131.12 


16 


7.125 


126.00 X 114.00 


0.015 


0.030 
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Table A. 3: Information for the three layers of proportional tubes in Station 4 as used 
in taking data for this analysis. Layers are listed in ascending distance from the 
origin. All lengths are in inches. 





Number 


Cell 


Aperture 




Operating 


Bank 






z 




of Tubes 


{x X y) 


{x X y) 




Voltage 


PT-Yl 


120 


1.00 X 1.00 


117.0 X 120.0 


2041.75 


+2500 


PT-X 


135 


1.00 X 1.00 


135.4 X 121.5 


2135.88 


+2500 


PT-Y2 


143 


1.00 X 1.00 


141.5 X 143.0 


2200.75 


+2500 
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Figure A.l: Block diagram of the E866 trigger system for the left hand side. Taken 
from JE] 
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Appendix B 
Analysis 



Table B.l: Estimated fraction of events in the target analysis of Run 2751 that 
originated in the dump. Dump IC3SB is the fraction of the luminosity that survived 
the target and interacted in the dump. Totals are for the three targets only. 



Target 


Total 


Dump 


Estimated 


Number 


Fraction 




IC3SB 


IC3SB 


Dump 


Events 


Dump 








Events 




Events 





40661 


40661 


143 ± 12 


143 


1.000 


1 


53013 


48708 


170 ± 15 






2 


52872 


46663 


163 ± 14 






3 


27225 


22983 


80 ± 7 






Total 


133110 


118324 


413 ± 22 


6350 ± 80 


0.065 ± 0.004 
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Table B.2: Spectrometer calibrations for data taken with the magnetic fields in SM12 
and SM3 parallel. All target analyses used Ztgt — —24.0 and Zg^at — 200.0 inches. 
All dump analyses used Zxg^ — 85.1 and Zgcat — 235.0 inches. 



Run 




i VV Ej Ej 


A U 


1 U 




CI' 

(rad) 












(rad) 


07/1 Q 


iarget 


u.yyo 


U.U4U 


-U.OZO 


U.UUUiO 


-U.UUUzo 




1 111 1 > "I T~\ 

VJ Llllip 


u.yyo 


u.u^u 


-u.ouu 


n nnnns 

U.UUUUo 


-U.UUUiO 


07/1 o 


Target 


n on 
U.yyt) 


U.Ut)i 


-U.OOO 


U.UUUiU 


-U.UUUoO 




1 In TTIT^ 

\-J Uilip 


u. yy 1 


u.uuo 




n nnnn^ 

U.UUUUO 


n nnnnn 

u.uuuuu 


2750 


Target 


0.990 


0.070 


-0.540 


0.00010 


-0.00035 




Dump 


0.992 


0.081 


-0.572 


0.00008 


0.00001 


2751 


Target 


0.990 


0.007 


-0.563 


0.00005 


-0.00035 




Dump 


0.992 


0.024 


-0.585 


0.00005 


0.00004 


2752 


Target 


0.990 


0.032 


-0.553 


0.00005 


-0.00035 




Dump 


0.992 


0.051 


-0.575 


0.00000 


0.00004 


2753 


Target 


0.990 


0.050 


-0.555 


0.00005 


-0.00035 




Dump 


0.992 


0.056 


-0.578 


0.00000 


0.00004 


2754 


Target 


0.988 


0.035 


-0.549 


0.00005 


-0.00035 




Dump 


0.990 


0.054 


-0.583 


0.00000 


0.00000 


2755 


Target 


0.986 


0.040 


-0.569 


0.00005 


-0.00035 




Dump 


0.990 


0.049 


-0.593 


0.00005 


-0.00015 


2756 


Target 


0.986 


0.049 


-0.545 


0.00001 


-0.00035 




Dump 


0.990 


0.072 


-0.581 


0.00000 


-0.00005 


2757 


Target 


0.986 


0.045 


-0.547 


0.00001 


-0.00035 




Dump 


0.990 


0.066 


-0.581 


-0.00002 


-0.00008 
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Appendix C 

Prom The Minimizations 

C.l Sub-Total And Total 
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Table C.l: The eight sub-total, total, number of degrees of freedom and reduced x^, 
Xpdf for all minimizations to the copper target data. Top line gives the individual 
spectra sub-totals, and the second line gives the sub-totals by charge and production 
region. 



Copper Data 


Function 


n 


2 2 2 2 
Xt+ X.T- ^1)+ 

X+ X- Xt Xd 


^c?/ Xldf 


Ai exp {—Bpt) 


NA 


18.43 21.68 16.90 24.41 
35.33 46.09 40.11 41.31 


81.42 61 1.33 




FF 


12.28 18.40 10.65 12.25 
22.94 31.16 31.19 22.91 


54.09 59 0.92 


4 


39.72 48.38 133.36 83.95 
173.08 132.32 88.10 217.31 


305.41 61 5.01 


[Pt + a m J 


5 


15.04 23.89 35.44 25.21 
50.48 49.10 38.93 60.65 


99.58 61 1.63 


6 


11.98 18.96 12.80 12.33 
24.78 31.29 30.94 25.14 


56.08 61 0.92 


A {^-PT / Pbeam)"' 
^2 / 2 1 2\" 

[Pj, + am^) 


4 


13.71 18.31 11.59 14.87 
25.30 33.18 32.03 26.46 


58.48 59 0.99 


5 


12.32 19.14 10.44 12.70 
22.76 31.85 31.47 23.15 


54.61 59 0.93 


6 


12.52 18.94 10.45 12.32 
22.96 31.26 31.46 22.76 


54.23 59 0.92 
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Table C.2: The eight sub-total, total, number of degrees of freedom and reduced x^, 
Xpdf for all minimizations to the beryllium target data. Top line gives the individual 
spectra sub-totals and the second line gives the sub-totals by charge and production 
region. 



BeryUium Data 


Function 


n 


2 2 2 2 

Xt+ Xt- Xd+ Xd~ 
Xl X- Xt Xl 


X^ ndf xldf 


Ai exp {—Bpt) 


NA 


9.55 7.44 12.71 33.18 
22.26 40.62 16.99 45.89 


62.88 63 1.00 




FF 


7.58 6.69 12.94 32.93 
20.52 39.62 14.27 45.87 


60.14 61 0.95 


4 


65.08 56.19 13.65 32.40 
78.73 88.60 121.27 46.06 


167.32 63 2.66 




5 


17.59 20.76 12.34 33.73 
29.94 54.40 38.36 46.08 


84.44 63 1.34 


6 


7.76 8.54 12.76 33.12 
20.52 41.66 16.30 45.88 


62.19 63 0.99 


4 (1 -Pr /P6eam)'" 
^2 / 2 , 2\" 


4 


7.28 6.15 12.91 32.96 
20.18 39.12 13.43 45.87 


59.30 61 0.97 


5 


6.69 5.88 12.91 32.96 
19.60 38.84 12.57 45.87 


58.44 61 0.96 


6 


6.59 5.79 12.90 32.97 
19.49 38.76 12.38 45.87 


58.25 61 0.95 
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